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Operator-sum or Kraus representations for single-mode Bosonic Gaussian channels are developed, 
and several of their consequences explored. The fact that the two-mode metaplectic operators 
acting as unitary purification of these channels do not, in their canonical form, mix the position and 
momentum variables is exploited to present a procedure which applies uniformly to all families in 
the Holevo classification. In this procedure the Kraus operators of every quantum-limited Gaussian 
channel can be simply read off from the matrix elements of a corresponding metaplectic operator. 
Kraus operators are employed to bring out, in the Fock basis, the manner in which the antilinear, 
unphysical matrix transposition map when accompanied by injection of a threshold classical noise 
becomes a physical channel, denoted T>(n) in the Holevo classification. The matrix transposition 
channels T>(k), T>(k~ 1 ) turn out to be a dual pair in the sense that their Kraus operators are related 
by the adjoint operation. The amplifier channel with amplification factor k and the beamsplitter 
channel with attenuation factor turn out to be mutually dual in the same sense. The action 
of the quantum-limited attenuator and amplifier channels as simply scaling maps on suitable quasi- 
probabilities in phase space is examined in the Kraus picture. Consideration of cumulants is used to 
examine the issue of fixed points. The semigroup property of the amplifier and attenuator families 
leads in both cases to a Zeno-like effect arising as a consequence of interrupted evolution. In the cases 
of entanglement-breaking channels a description in terms of rank one Kraus operators is shown to 
emerge quite simply. In contradistinction, it is shown that there is not even one finite rank operator 
in the entire linear span of Kraus operators of the quantum-limited amplifier or attenuator families, 
an assertion far stronger than the statement that these are not entanglement breaking channels. A 
characterization of extremality in terms of Kraus operators, originally due to Choi, is employed to 
show that all quantum-limited Gaussian channels are extremal. The fact that every noisy Gaussian 
channel can be realised as product of a pair of quantum-limited channels is used to construct a 
discrete set of linearly independent Kraus operators for noisy Gaussian channels, including the 
classical noise channel, and these Kraus operators have a particularly simple structure. 

PACS numbers: 03.67-a, 03.67.Mn, 03.67.Hk, 42. 50. Ex, 42.50.-p 



Gaussian states and Gaussian channels play a major role in quantum information processes, and this is primarily 
due to their potential experimental realization within current technologies [H-Q- Indeed, the basic protocols of 
quantum information processing including teleportation and dense coding have been implemented in the quantum 
optical domain @, @. The feasibility of processing information using Gaussian channels was originally explored in 
[l|, [7[ . More recently, the problem of evaluating the classical capacity of Gaussian channels was addressed in 



and the quantum capacities in [llMlq . In particular, the classical capacity of the attenuator channel was evaluated in 
[l(| , and the quantum capacity of a class of channels was studied in [l2| • A systematic study of the structure of the 
family of all Gaussian channels has been carried out in (l7r{2l1 |; single- mode Gaussian channels have been classified 




Phase space description in terms of quasiprobabilities or, equivalently, in terms of the associated characteristic 
functions underlies the very notions of Gaussian states and Gaussian channels. Gaussian states are those with 
Gaussian characteristic functions, and Gaussian channels are those trace-preserving CP maps which image every 
input Gaussian state into a Gaussian state at the output. 




I. INTRODUCTION 
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Gaussian states are fully specified by their first and second moments. Since the first moments play no significant 
role in our study, we may assume that they vanish (this can indeed be ensured using the unitary Weyl-Heisenberg 
displacement operators), so that a Gaussian state for our purpose is fully described by its covariance matrix [22M25j | . 
The symplectic group of real linear canonical transformations (acting through its unitary metaplectic representation) 
and the Weyl-Heisenberg group of phase space translations are the only unitary evolutions which preserve Gaussianity, 
and these groups are generated by hermitian Hamiltonians which are respectively quadratic and linear in the creation 
and annihilation operators (23 - l24| . This suggests that more general Gaussian channels on system A may be realized 
as Gaussianity preserving unitaries on a suitably enlarged system : 

PA -> p' A = Trs (u A b {pa ® Pb)u\ b ) . (f.l) 

Here pb is a Gaussian state of the ancilla B, and Uab is a linear canonical transformation on the enlarged composite 
system consisting of the system of interest A and the ancilla B. That all Gaussian channels can indeed be realized in 
this manner has been shown by the work of Holevo and coauthors [13, EH, US HJ . 

It is clear that the most general trace-preserving linear map ft which takes Gaussian characteristic functions to 
Gaussian, taking states with vanishing first moments to ones with vanishing first moments, are necessarily of the form 
^ : X(£) ~~ ► X (£) = x{X£,) exp[— i£ T Y£], where X, Y are real matrices with Y — Y T > 0. And X, Y need to obey 
an appropriate matrix inequality to ensure that the trace-preserving map Q is completely positive (Til fl2l I2H |27| . 
For a given X, the minimal Y, say Yq, meeting this inequality represents the threshold Gaussian noise that needs to 
be added to x(A£) to make atonement for the failure of A to be a symplectic matrix, and thus rendering the map 
completely positive; if X happens to be a symplectic matrix, then the corresponding minimal Yq = 0. 

Now, given a Gaussian channel f2 we can construct, 'quite cheaply', an entire family of Gaussian channels by simply 
preceding and following fi with unitary (symplectic) Gaussian channels U{Si),U{S 2 ) corresponding respectively to 
symplectic matrices Si, S 2 . Therefore in classifying Gaussian channels it is sufficient to classify these orbits or double 
cosets and, further, we may identify each orbit with the 'simplest' looking representative element of that orbit (the 
canonical form). Since 

U{Si)SW{S 2 ) : x(0 -^ X {S2XSi0exp[-±fSlYSiZ], (1.2) 

the task actually reduces to enumeration of the orbits of {X, Y) under the transformation (A, Y) —> (A , Y ) — 
{S 2 XS U SfYSi). 

One final point before turning to the special case of single-mode Gaussian channels. The injection of an arbitrary 
amount of classical (Gaussian) noise into the state is obviously a Gaussian channel : %(£) — > x{Q ex P[ — §£ T £L a > 0. 
It is called the classical noise channel. Now, given a Gaussian channel we may follow it up with a classical noise channel 
to obtain another Gaussian channel. A Gaussian channel will be said to be quantum-limited if it cannot be realized 
as another Gaussian channel followed by a classical noise channel. Conversely, the most general Gaussian channel is a 
quantum-limited Gaussian channel followed by a classical noise channel, and it follows that quantum-limited channels 
are the primary objects which need to be classified into orbits. 

In the single-mode case where A, Y are 2x2 matrices, Si, S 2 £ Sp{2, R) can be so chosen that A equals a multiple 
of identity, a multiple of (73, or (1 ~t~ 0*3 )/2 while Y equals a multiple of identity or (1 + 0-3) /2. Thus the canonical 
form of a Gaussian channel A, Y is fully determined by the rank and determinant of A, Y and we have the following 
classification of quantum-limited bosonic Gaussian channels (l7l . EH 



P(k;0) 
Ci(k;0) 
C 2 (k;0) 
Ai{0) 
A 2 {0) 
B 2 {0) 
Bi(0) 



X = -n<7 3 , Y = (1 + k 2 )1, k > 0; 

X = k1, Y = (1 - k 2 )1, < k < 1; 

X = k1, Y q = {k 2 - 1)1, k > 1; 

A = 0, Y = l: (1.3) 

A = (1 + cr 3 )/2, Y = l; 

A = 1, Y = 0; 

A = 1, Y = 0. 



It may be noted that the quantum-limited end of both the Bi and B 2 families is the trivial identity channel. 

By following the above listed quantum-limited channels by injection of classical noise of magnitude a we get 

respectively T>{n;a), Ci{n;a), C 2 {k;o), Ai{a), A 2 {a), and B 2 {a); the last case £>i(a) is special in that it is obtained 

from Z?i(0) by injection of noise into just one quadrature: x{Q ~> x{Q ex P[ — f£ T (l + "^Ki- 
lt is clear in the case of T>{k; 0) that A = — ko~j, corresponds to (scaled) phase conjugation or matrix transposition of 

the density operator. And the phase conjugation is the most famous among positive maps which are not CP (3214341 ]: 
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it is the injection of additional classical noise of magnitude (not less than) 1 + k 2 , represented by Yq, that mends it 
into a CP map. 

The reason for the special emphasis on quantum-limited channels in our enumeration of the Holevo classification 
is this : every noisy Gaussian channel [except Bi(a)] can be realized, as we shall see later, as the composite of a pair 
of quantum-limited channels. And this fact proves to be of much value to the study presented in this paper. In 
the original classification of Holevo [13, E3 the families C\ and C2, which correspond respectively to the attenuator 
(beamsplitter) and the amplifier (two-mode squeezing) channels, together constituted a single family C. From the 
point of view of the present study, however, these two families turn out to be qualitatively different from one another, 
hence we prefer to keep them as two distinct families. 

It is well known that every trace-preserving completely positive map has an operator-sum representation of the 
form 



often called Kraus representation [28|. It may be noted, however, that this representation appears as Theorem 4 
of a much earlier work of Sudarshan et al (29[. It has been presented also by Choi [3(J, apparently independently. 
Mathematicians seem to view it as a direct and immediate consequence of the dilation theorem of Stinespring (3"lj . In 
this paper we develop and present a systematic analysis of the operator-sum representation for single-mode bosonic 
Gaussian channels. 

Knowledge of Kraus representation of a channel could prove useful for several purposes. For instance, since the set 
of channels for a given system is convex, it is of interest to know its extremals. And a theorem of Choi [3(j gives 
a necessary and sufficient test for extremality of a channel in terms of Kraus operators. It is known that a channel 
is entanglement breaking if and only if it can be described in terms of a set of rank one Kraus operators [35l - f37l | . 



Further, the work of [38( and [39|, |40| relate error correctability to the structure of the Kraus operators of a channel. 
Finally, there has been considerable recent interest in contrasting the Gaussian with nonGaussian states in respect 
of robustness and degradation of bipartite entanglement under one-sided and two-sided action by Gaussian channels 
[4T| - [43l | , and it is likely that Kraus representation could throw light on this problem. 

The content of this paper is organised as follows. Section II presents a general scheme for computation of Kraus 
operators, and this scheme applies uniformly to all quantum-limited Gaussian channels. This scheme takes particular 
advantage of the fact that the symplectic two-mode transformation which realizes the channel in the sense of 
does not couple, in the Holevo canonical form, the position variables with the momentum variables. With the ancilla 
mode assumed to be in its vacuum state initially, it turns out that the Kraus operators for each channel can be simply 
read off from the matrix elements of the appropriate two-mode metaplectic operator. 

This scheme is applied in Sections III to VI to detail the Kraus operators of respectively the T>, C\, C2, and A2 
families of quantum-limited channels. The Kraus operators in every case is found to have an extremely simple-looking 
sparse structure. In each case we ask if the channel has any fixed points (invariant states), and if there are sufficient 
number of rank one operators in the linear span of the Kraus operators, a question which is at the very root of the 
entanglement breaking property of the channel. 

In the case of the phase conjugation or (matrix) transposition family T>(n; 0) in Section III we explore how the 
threshold noise of magnitude 1 + k 2 renders the antilinear phase conjugation into a linear map and channel. We 
bring out a well-defined sense in which the channels T>(k; 0) and 2?(k -1 ; 0) are dual to one another. The case T>(1; 0) 
is self-dual and hence doubly stochastic, but it turns out that it is not a random unitary channel, a fact which is 
of relevance to the possibility or otherwise of extending the classical Birkhoff theorem to the quantum domain [52| . 
We examine in the Kraus picture the manner in which quantum- limited attenuator and amplifier families Ci, C2 act 
as simple scaling maps respectively on the diagonal weight function and the Husimi Q-distribution. Comparing the 
Kraus operators of the C\ family with those of the C2 family, we show in Section V that these two families are dual 
to one another. The intersection between C\, C2 consists of just the identity channel, the only self-dual or doubly 
stochastic channel in the union of Ci, €2- The manner in which the semigroup structures of the C\ and C2 families 
are reflected in their respective sets of Kraus operators is brought out in Sections IV and V, and this enables us to 
point to a Zeno-like effect [55| . in both cases, arising as consequence of interrupted evolution. Finally, even though the 
single-quadrature classical noise channels Si (a), a 7^ [£>i(0) is the identity channel] are not quantum- limited, we deal 
with them briefly in Section VII just to bring out the fact that this case too is obedient to the general computational 
scheme presented in Section II. 

In Section VIII wherein we use Choi's theorem [3(| to study if there are any extremals among Gaussian channels, 
we show that all quantum limited channels, and these alone, are extremal. That our concern upto this stage of the 
presentation is (almost) exclusively with the quantum-limited case gets justified by our demonstrations in Section IX 
that every noisy Gaussian channel [except £>i(a)] can be realized as the composite of a pair of quantum-limited ones. 
This demonstration leads, in particular, to an operator-sum representation for all noisy channels, including the all 




(1.4) 



a 



a 
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important classical noise channels 62(a) [but excluding Bi(a), and only this case], in terms of a linearly independent 
discrete set of Kraus operators having very simple sparse structure. 

The composition of pairs of quantum-limited channels studied in Section IX, and conveniently summarized in Table 
U there, assumes that both the constituent channels are simulatneously in their respective canonical forms. When 
this assumption is removed, the situation with the composition process gets much richer. The general case is fully 
classified and presented in Table [TT| of the Appendix. 

The final Section X contains a brief summary of the principal results and also some additional remarks. 



II. KRAUS REPRESENTATION: SOME GENERAL CONSIDERATIONS 

Given density operator p^ describing the state of a single-mode radiation field, the action of a quantum-limited 
Gaussian channel takes it to [l?! [HI 

= Tr b {U^ (p^ <g> |0) b6 (0|) t/ (ab)t ). (2.1) 

Here 1 0) t, is the vacuum state of the ancilla mode b, and f/' ab ' is the unitary operator corresponding to a suitable 
two-mode linear canonical transformation. It is convenient to perform the partial trace in the Fock basis of mode b. 
We have 

p'M = J2b{t\U (ab) (P (a) ® |0) 6b (0|)(7 (a6)t \£) b 

i 

= J2bm {ab) \o)bP {a) b (o\u^ b ^\e) b . (2.2) 

t 

Clearly, b (£|£/( ab >|0) b is an operator acting on the Hilbert space of mode a. The last expression thus leads us to the 
Kraus representation of the channel [28| : 

P^P (a) =J2 w fP (a)w L W e = b (£\U^\0) b . (2.3) 

e 

It follows that once the Fock basis matrix elements of U^ ab ' are known, the Kraus operators Wi can be easily read 
off. Let {mim 2 \U^ ahS> \nin 2 ) = C™ 1 ™™ 2 be the matrix elements of [/( ab ) in the two-mode Fock basis. Since the ancilla 
mode b is assumed to be in the vacuum state, the W/s are obtained by setting n 2 = and m 2 = £ : 

00 

w *= E <CoVi><ni|- (2.4) 

Tl\ ,7711=0 

Now, in evaluating C™^ 12 it proves useful to employ a resolution of identity in the position basis [44[ : 

C™^ = ( mim2 \uW\n in2 ) 

poo 

dx\dx 2 {vn\m 2 \x\x 2 ) (x\X2\U { ' ah ^ \ri\ni). (2-5) 

Under conjugation by \J^ ah ^ the quadrature variables qj,Pj (j = 1,2) undergo a linear canonical transformation 
S € Sp(4,R), of which U^ ah ^ is the (metaplectic) unitary representation J24[. Let us assume that this canonical 
transformation does not mix the position variables with the momentum variables. That is, 

Qi) _+ u wi(Qi) u w = (li\ =M ( qi 

& J \12) \q 2 ) \q2 

where M is a real non-singular 2x2 matrix. This assumption that our S E Sp(4, R) has the direct sum structure 



(M ) r L ), (2.6) 
P2J \P 2 J \P 2 J \P2j 



5 



S = M © (M ) will prove to be of much value in our analysis. We have 



nin 2 



dxidx2{mim2\xiX2){xix 2 \U { - ab ' > \n 1 n2) 
dx 1 dx2(m 1 m2\xiX2)ipn 1 (x' 1 )ip n2 (x , 2) 
dx 1 dx 2 ip* ll (xi)ip^, 2 (x 2 )Vv (x[)ip n2 (x' 2 ), 



(2.7) 



where (x 1; x 2 ) is linearly related to (xi, X2) through M. These wavefunctions are the familiar Hcrmite functions, the 
Fock states in the position representation. The above integral may be evaluated using the generating function for 
Hermite polynomials [44| : 



-1/4 



n -l/4 gn 



e- x H n (x) 



1 



1pn(x) 



Inserting in Eq. (|2.7[) the generating function for each of the four wavefunctions we have 



exp[--[(x-zV2) 2 -z 2 } 



Qm, 1 Qni2 Qn^ Qn 2 



V^Wm^l dr,™ 1 dry™ 2 dz^ dz. 



zi,z 2 ,r) 1 ,ri2=0 



(2.8) 



(2.9) 



where 



F(z 1 ,z 2 ,r) 1 ,r) 2 ) = tt~ 



dxi dx2 exp 
+(xi - Zl V2) 2 + (x' 2 - Z2V2) 



-[(xi - r? x V2) 2 + (x 2 - 772 V2) 2 



22222 



}■ 



(2-10) 



The Gaussian integration over the variables X\ and X2 can be easily carried out to obtain F(z\, z 2 , 171, ^2)1 an d from 
F(zi, z 2 , 771, 772) we may readily obtain C™ 1 ^ 2 , and hence the Kraus operators. This is the general scheme we will 
employ in what follows to obtain Kraus representation for quantum-limited Gaussian channels of the various families. 



III. PHASE CONJUGATION OR TRANSPOSITION CHANNEL V(k), k>0 

We now use the above scheme to evaluate a set of Kraus operators representing the phase conjugation channel. 
The metaplcctic unitary operator U^ ab ^ appropriate for this case induces on the quadratur e op erators of the bipartite 
phase space a linear canonical transformation corresponding to the following S € Sv(4l, R) [17| : 



'sinh/i cosh /1 

— sinh fi cosh \i 

cosh fi sinh fj, 

cosh /1 — smhfi/ 

Written in detail, the phase space variables undergo, under the action of this channel, the transformation 



(3.1) 



w \q 2 J V?2 



I) - (2) - <W - ,,T (S 

sinh fi cosh fi 



^ I cosh /i — sinhjuy ' (3-2) 

It is seen that the above 5 is indeed of the form S = M© (M _1 ) T € 5p(4, R), and does not mix the position variables 
with the momentum variables, and so our general scheme above readily applies. 
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It is clear from the structure of S that the parameter /i is related to k in T>(k) through n — — sinh/i > 0, so that 
cosh /i = V ' n 2 + 1. Thus (|2.10[) translates, for the present case, to the following expression: 

F(zi,z 2 ,T] 1 ,r] 2 ) = 7r _1 J dxi dx 2 exp j - - ?/iV2) 2 + (x 2 - 772 V2) 2 

+(-KXt + \J\ + k 2 x 2 - ziV2) 2 + {s/l + k 2 xi - ra 2 - z 2 \/2) 2 
-^-% 2 -^-^]|- (3-3) 

Performing the Gaussian integrals in x\ and x 2 we obtain 

F{zi,z 2 ,rji,rj 2 ) = (yj 1 + k 2 ) -1 exp|(i/l + k -2 ) -1 (771772 - ^1^2) 

+ (-v/l + K 2 )- 1 (77^2 + 772Z!)}. (3.4) 

To obtain the matrix elements C™^™ 2 we need to carry out the procedure indicated in Eq. (|2.9[) . This may be 
done in two steps. We begin by rewriting the function F{z\,z 2 , 771,772) as 

F(zi, z 2 , 771, 772) = (V« 2 + l) _1 cxp{z 2 [( v / l + K 2 )- 1 77i - (V 1 + k- 2 )- 1 z{\ 

+772 KV 7 ! + k- 2 ) _1 t7i + (Vl + k 2 ) _1 zi]} . (3.5) 
Performing the 2:2 and 7/2 differentiations respectively 712 and 7712 times on F{z\, z 2 , 771, 772), we obtain 



[(Vl + ^rV - (Vl + k- 2 )- 1 ^]" 2 ^^! + k- 2 ) _ V + (Vl + ^r^iP 2 ^ = GF. (3.6) 

The remaining differentiations can be carried out using the Leibniz rule. Since we finally set 21,22,771,772 = 0, and 
since F(0) = 1, the only terms that could possibly survive are necessarily of the form 



^m-r^rKVi+^TV - (Vi + ^r^ipKVi + ^rv + (VT+^r i z 1 ] m \ (3.7) 



To evaluate the above expression we set x = (V k 2 + 1) 1 771 — (vT+ k~ 2 ) and y = (VT + k~ 2 ) 1 
(Vl + k 2 )~ 1 zi, and compute 



'/1 



[(v/1 + k 2 )- 1 ^ + (VI + K-^dy}^ [-(Vl + ^)- 1 d x + (Vl + K 2 )" 1 ^]" 1 x" 2 j7 m2 U !y=0 . (3.8) 
Straight forward algebra leads, in view of Eq. (|2.9[) . to 



erf = (^ 1 + K2 ) _ y" y»ic, mi a (-Vi + K- 2 )~ (mi+j ' _r) (Vi + K 2 )" (ni ^ +r) 

1 2 V77l!772!mi!TO 2 ! ^ 

x (-i)™i- r ra 2 !m 2 !(5 n2ir+ j 5 TO2jni _ i+TO1 _ r . (3.9) 

The Kraus operators We, denoted X^(re) in this case, are obtained from these matrix elements by setting 712 = and 
7n 2 = t- Since ?i 2 = =>■ r, j = 0, we have, 



2 \— 1 ™i m i 



t £ (k) = (VT+^r 1 V (v r +^)'" 1 (^ 1 + ^ 2 )' T " 1 v^ ^ ni+mi( „ iri [mi)<ni|- (310) 

We set 7ii + mi = £ and denote 7ii = n, leading to 



T/(k) = (Vl + « 2 r : X)(V1 + « 2 )""(Vl + K- 2 )- ( '- n) V^K - n)(n|, £ = 0,1,2,- •• (3.11) 



71=0 



as our final form for the Kraus operators of the phase conjugation channel. We note that the T^(k)'s are real and 
manifestly trace-orthogonal: tr(Tg(«;)TT^/ («;)) = if £ 7^ £ . 
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I. The dual of V(k) 

As is well known (and also obvious), if a set of Kraus operators {We} describes the completely positive map 
$ : p — > p = ^2iWipWg, then the dual map $ : p — > p = ^gWi' pWi, described by the dual or adjoint set of 
operators {W ( , }, is also completely positive. It is clear that the dual map $ is unital or trace-preserving according as 
3" is trace-preserving or unital. 

For the present case of T>(k), it is readily verified that the Kraus operators {Tf(re)} presented in (|3.1ip meet 
^2gTe' (k)T^(k) = 1, consistent with the expected trace-preserving nature of p — > p = ^2 e Te(n) pTe^ (k). But the 
phase conjugation channel is not unital in general, for we have 

£}t<(k) = k" 2 1. (3.12) 

i 

We may say that it is 'almost unital' to emphasise the minimal nature of the failure : the unit element is taken by the 
channel into a scalar multiple of itself. However, the scale factor k~ 2 can not be transformed away by absorbing k _1 
into the Kraus operators, for the Kraus operators so modified would not then respect the trace-preserving property 
of the map. 

It is thus of interest to understand the nature of the unital channel described by the set of Kraus operators {T^ny}. 
We have 

l 

T e ( K y = (Vl + k 2 )- 1 + + K- 2 r {e ~ n) V^|n><* - n\ 

n=0 


= (v/IT^r 1 J2 (a/i+^T^'^a/i + K- 2 r n ' V 1 ^^ \l - n')(n'\ 




= (Vl + K 2 )- 1 J2(^ k2 + ir {e - n) {\Tl + k-^V^M - n)(n\ 

= k^T^k,- 1 ). (3.13) 

Thus the dual {Te(nY} differs from the original {X^(k)} in two elementary aspects. The multiplicative factor kT 1 is 
the same for all Kraus operators, independent of £. Thus the only significant difference is change in the argument of 
Ti, from k to We conclude that the 'dual' channel whose Kraus operators are k 2i(/c)t is the (trace-preserving) 
phase conjugation channel I?(k~ 1 ). We have thus proved 

Theorem 1 While the Kraus operators {Te(n)} describe T>(k), the 'dual' channel described by Kraus operators 
{kT((kY} is the trace-preserving phase conjugation channel 2?(k -1 ) with reciprocal scale parameter. 



II. Properties of the Kraus operators 

We now explore the properties of the Kraus operators of Eq. (|3.1ip , connecting their explicit action in the Fock 
basis to the expected transformation of the characteristic function. The question of its fixed points is studied through 
the action of the channel on the cumulants. The action of the channel is illustrated with simple examples and, finally, 
the entanglement breaking nature of the channel is made transparent by obtaining a set of rank one Kraus operators 
describing the channel. 

The expected or defining action of the phase conjugation channel on the characteristic function is fl7j : 

Xw(0 -> Xw(0 - Xw(-«0 exp[-(l + K 2 )\£\ 2 /2}. (3.14) 

It is of interest to understand how the 'antilinear' phase conjugation (£ — > £*) action of this channel on the character- 
istic function emerges from the linear action of the Kraus operators. To this end, it is sufficient to establish such an 
action on the 'characteristic function' corresponding to the operators |n)(m|, for arbitrary pairs of integers n, m > 0. 
The 'characteristic function' of |h)(?ti| is given by [45| 

Xw\ n ){m\(0 = ( m \ D (0\ n ) 

" ro (|£| 2 ) exp[-|£| 2 /2] for n > m, 
" n (|C| 2 ) cxp[-|e| 2 /2 for n<m. (3.15) 
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Assuming n > m, the action of the phase conjugation channel on the operator |n)(m| is 

oo 

^T e (K)\n)(m\Tj(K) 



= (1 + k 2 )- 1 + «2)-( n+m i (Vl + K-2)-( 2 '-"- m > y/'C n e C m \£ -n)(£- m\. 

Denoting n = m + d and ^ — n = A, we have 



(3.16) 



Y,T t {n)\m + 5){m\Tl{K) = (1 + n 2 )- 1 ^! + K 2)-(2m+«) ^ + K _ 2 )-a 



^ n /(to + 5)!to!A!(A + (5)! 



(3.17) 



The manner in which T>(k), matrix transposition accompanied by threshold Gaussian noise exp[— (1 + k 2 )|£| 2 /2], acts 
as a channel may now be appreciated. Every operator M can be written in the Kronecker delta basis as 
M = y] j £ Cjt\j){£\. The coefficient matrix C associated with |5)(3|, for instance, is Cj^ = SsjSes, with non-zero entry 
only at the lower-diagonal location (5, 3) marked £§> in the matrix below. 



/0 0x00000 
000x0000 

ooooeooo 

00000x00 
x0 

OOOOOOOx 
00000000 



V 



/ 



On transposition this entry moves to the upper-diagonal location (3, 5) marked ©, and the threshold noise then spreads 
it along the parallel upper diagonal (3 + r, 5 + r), — 3 < r < oo marked x . 

Let the Weyl-ordered characteristic function tr(D(£)|m + S)(m\) where £>(£) = exp[£a^ — £*a] is the displacement 
operator, be denoted Xw\ m +5}{m\(Q> an< ^ that °f tne output J27Lo Ti{K)\m + 6) (m\Tg(Ky be denoted Xwi m +$\ /„,](£)• 
Then we have from Eq. (|3.17p 



xw' lm+s)(m \(0 = (i + n 2 )- 1 (^/TT^r {2m+s \VT+^) 

(A + m + 5)!(l + K- 2 )-* 



-6 



E 



A ^ v/(m + <J)!m!A!(A + 5)! 
(1 + K 2 )- x e-\^ 2 ' 2 



■(X + S\D(0)\X) 



y/(m + S)lm\ 

OO 



( Vl + K 2 ) - (2 " i+5) ( Vl + K- 2 ) 



2^ — <5 



A=0 



(A + 5)! 



(3.18) 



where we used (|3.15|) . the Fock basis representation of the displacement operator. While no 'phase conjugation' is 
manifest as yet, we expect from Eq. p,14[) that the channel should take the characteristic function of \m + 8)(m\ to 



Xw\ m+S )(m\(0 = (m\D(-K£*)\m + 6) exp 
= (m + 5\D(nt;*)\m)* exp 



--(1 + « 2 )|CI= 



;(1 + « 2 )ICP 



(nO S L s m (K 2 \i\ 2 ) exp 



m + S\ 



(3.19) 



9 



Thus the problem reduces to one of establishing equality of X W \ m +S)( m \(0 in P- 18 l> and Xw\ m +5)( m \(0 in <I 3 - 19 P- 
That is, it remains to prove 



y^y(«0 i im(« 2 |^| 2 )ex P [-(l/2 + « 2 )|C| 
(l + K 2 )-^"^ 2 / 2 



£(i + «- 2 )- A (yT+^)- (2m+5) (yT+^)- 



sfijn + S)lm\ 



(X + m + S )\ &STSn 

{x + sy. 



A=0 

eLim 2 ), (3-20) 



for all m, S > [the case of |m)(m + 5\ can be handled similarly]. 

Since the associated Laguerre functions form a complete orthonormal set, we may expand the LHS of Eq. (|3.20[) 

in the Laguerre basis. That is, we multiply both sides of Eq. (|3.20|) by (C*) S Lf(\£,\ 2 ) e~^~ and evaluate the overlap 
integrals. We use the following two standard results \ Ji) orthogonality relation among Laguerres, and (ii) the overlap 
between a Laguerre and a scaled Laguerre function |46[ : 

e^« |2 |?| 2 ^(IC| 2 )^(l^| 2 )rf|e| 2 =^^^,™. 



~ m * \e S Li(v 2 \Z\ 2 )L S M\ 2 )d\tf 



(m + e + Sy. (t-rf) m (t- 1) 



x F 



-m, —m — £ — 6, 



t(t - rf - 1 



(t-i)(t-n>) 



(3-21) 



Here F[-} is the hypergeometric function. In our case t = r/ 2 + 1, which implies that the last argument of F[-] in Eq. 
(|3.21[) is zero, and thereby F[-] = 1. Performing the overlap integrals, we obtain for the left and right hand sides of 
(13^201) 

LHg _ (m + £ + 6y. k 2 *+ s 



e\y/(m + S)\m\ (1 + ' 

RHS = ( w + £ + <5)! (VTT^)- (2+2m+ V (VT+^y^ +5 \ (3.22) 
<V(m + 6)!m! 

These two expressions obviously equal one another for all I. We have thus established Eq. (|3.20p . and the fact that the 
Kraus operators indeed effect the 'completely positive phase conjugation' operation, transforming the characteristic 
function as expected in ()3.14p . 

/ 1 1 1 2 

Theorem 2 The scaled phase conjugation transformation Xw(0 ~ * X w(£) ~ Xw(~ K £,*) exp[— (1 + k2 )2~] * s > ^ n 

view of the threshold noise exp[— (1 + k 2 )|£| 2 /2] a completely positive map, and is implemented linearly by the Kraus 

operators {Tz(n)} in Eq. (|3.11[) . 

The phase conjugation channel has an interesting property in respect of classicality/nonclassicality of the output 
states. We may say a channel is nonclassicality breaking if the output of the channel is classical for every input state. 
That is, if the normal-ordered characteristic function Xn(0 °f the output, related to the Weyl-ordered characteristic 
function Xw(Q °f (13.141) t hrough Xn(0 = Xw(Q ex P[l£| 2 /2], is such that its Fourier transform, called the diagonal 
'weight' function <fi(a) [471 ] . is a genuine probability density. 

Now, Eq. (|3.14j) written in terms of the normal-ordered characteristic function reads 

Xn(0 -► Xn(0 = Xw(-/<*)exp[- K 2 ri72] 

= Xa(-k?), (3-23) 

where xa(C) = Xn(0 ex P[ — 1£| 2 ] is th e antinormal-ordered characteristic function corresponding to the Q or Husimi 
distribution. 

Under Fourier transformation this important relation (|3.23[) , namely Xat(0 = Xa{— «^*) ; reads that the output 
diagonal weight function <fi (a) evaluated at a equals the input Q(a) evaluated at k a*. Thus <fi (a) is a genuine 
probability density for every input state, and we have 

V[ K ) : in (a) -> out (a) = K - 2 Q in ( K "V). (3.24) 
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Since the Q-distribution of a density operator is given by Q(a) = (a\p\a), it is a genuine probability distribution for 
all states including nonclassical states. We have thus proved 

Theorem 3 The phase conjugation channel is a nonclassicality breaking channel. 



III. Fixed points 

We now study the fixed points of the phase conjugation channel through consideration of cumulants. The char- 
acteristic function of the s-ordered quasiprobabilities differ from one and another just by a Gaussian factor, and it 
follows that the cumulants of order > 2 of the s-ordered quasi-probability are independent of the ordering parameter 
s [3 [4i|, and are thus intrinsic to the state. Hence it is sufficient to work with a particular choice of s. We work 
with s = 0, the case of symmetric or Weyl ordering. 

Given a symmetric ordered characteristic function xw(£)> the corresponding cumulant generating function is defined 

as 

r(0 = iog[xw(01- (3-25) 

With £ = £i + i^2 , the cumulants are defined through 

Tmi.m, = g {i ^ mi d(^ 2 )m 2 r (£)le=0, mi, 012 = 0, 1,2,." (3.26) 
From Eq. (|3.14[) we know that under the action of the phase conjugation channel 

V(k) : r(0 -> r'(0 = log [Xw(0] = log bcw(-^*)) ~ ^l 2 (i + (3.27) 



are 



Since the additional term on the right hand side is quadratic in £, the cumulants 7 mo of Xvy(£) °f order mi, 7^ 2 

_9_\ mi /_a 
0(w we 



7™ = ( tJttV 11 (-J^) 



5=0 



Denoting — k£* = t we have 



■mo 



7„, TO - (-irer* (s^j) (s^jj Kh^wi 

= (-lrW""'!.,,.,, (3-29) 

Now, for a state to be invariant all its cumulants need to remain invariant. By p.29|) . none of the cumulants of 
order mi, > 2 are preserved for k / 1. Indeed, under repeated use of the channel the higher order cumulants 
monotonically increase or decrease depending on whether k is > 1 or < 1. In the case k = 1, the cumulants with mi 
or ni2 = 2 are not preserved because of the last additional term on the right hand side of (|3.27[) . showing that no 
non-Gaussian state can be a fixed point of T>(k). 

We are therefore left with the case k < 1 to consider. It is clear that in this case all cumulants of order ^ 2 die out 
under repeated use of the channel, and any initial state is driven towards a fixed Gaussian (thermal state). A similar 
situation was discussed in |5(j, where linear devices were used to drive non-Gaussian pure states to Gaussian states. 

We are thus led to look for fixed points among Gaussian states. The additional last term in Eq. Q3.27P is proportional 
to |£| 2 . Since there is no cross term involving the real and imaginary parts of £, it is sufficient to look for fixed point 
among the thermal states, given by 



Pth(ao) = " , ( 1 1 ) l n )H 



(l-a05>»R (3.30) 



where x = (ao — l)(ao + 1) 1 , an d the average photon number tr (p^a'a) = (do — l)/2. By Eq. (|3.17j) the output of 
the channel is 



00 00 



p = {l + K 2 )- l {l-x)^^ n+j C j (l + K 2 )- j {l + K- 2 )- n x i \n){n\. (3.31) 

n=0 j=0 
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No. of applications 

FIG. 1: Showing the variation of the output thermal parameter a under repeated application of the quantum-limited phase 
conjugation channel T>{n) with k = 0.8 for various input thermal parameters. The stable thermal parameter is ao = 4.56. 



With the use of the identity 



oo 

g^QMl + *>)-V = (1 _ (1 + w2) _ la;)n+1 , (3-32) 



the double summation in Eq. (|3.31[) reduces to a single sum, and we have 

1 — x — ' 



71=0 



Pth(oo) -> Pth(Oo) = i + K 2 X _ x J2 f 1 + h K 2_ x ) l n )H ( 3 ' 33 ) 



This is a thermal state, and we see that the effect of the channel is to change the thermal parameter x as follows : 

k 2 

x^x = 5 . (3.34) 

1 -\- K — X 

Thus the fixed points are x = k 2 and x = 1. The first, x = k 2 , corresponds to a finite temperature state (k 2 < 1), 
to which all other states are attracted under repeated use of the channel. The second one corresponds to infinite 
temperature, and is uninteresting for this reason and for the fact that no state is attracted towards it. 
In terms of the parameter ao, the channel action reads 

V(k) : pth(ao) -> Pth(a ), a Q ^ a = K 2 a Q + l + K 2 . (3.35) 

As was to be expected from (|3.34[) . the recursion relation in Eq. (|3.35[) is stable when k 2 < 1 and unstable for n 2 > 1. 
For a given k < 1, any input state is driven towards this 'stable' thermal state under repeated use of the channel, and 
the thermal parameter of this attracting fixed state is ao = (1 + k 2 )(1 — k 2 )^ 1 . This is illustrated in Fig. (JXJ) . We 
may summarize as follows 

Theorem 4 For the phase conjugation channel T>(k), k < 1, there is a unique thermal state /Oth(ao) with parameter 
oq = (1 + n 2 )(l — k 2 ) -1 which is left invariant by the channel. All other states are driven towards this thermal state 
under repeated use of the channel. Channels T>(k) for which k > 1 have no fixed points. 

As a simple illustration of the action of the phase conjugation channel, assume the input to be a Fock state. By 
Eq. ([5Tf]) . we have 



\n)(n\^yT £ (K)\n)(n\T e (K^ 



(1 + k 2 )- 1 e °n(l + k 2 )""(1 + K- 2 )-V- n) \l -n){£-n\. (3.36) 
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Setting I — n — j, we have 

j£T e (K)\n){n\Tt( K )1 = (1 + re 2 )" 1 E" + ^-(l + « 2 )""(1 + O" j b')01- (3-37) 

That is, a Fock state is taken to a convex sum of aZZ Fock states. As an immediate consequence we have : if a density 
operator p has (n|p|n) = for some n, then p cannot remain invariant under the action of the channel. This is true, 
in particular, of any state p which is in the support of a finite number of Fock states. 
As another simple example, consider the phase averaged coherent state given by 

P = e- X J2^\j)Ul (3-38) 
This mixed state has Poissonian photon number distribution [PND]. Under the action of the channel, we have 

OO , j oo 

p P = e-\l + re 2 )" 1 J2 7, E n+3C ^ 1 + "V(l + K~T»<4 (3-39) 

We cannot solve consistently for parameters re and A such that the output PND is also Poissonian. In view of 
Theorem [3J the output is a classical state, and hence is necessarily super-Poissonian (Hlj . To summarize, under the 
action of the channel T>(k) a thermal PND is taken to a thermal PND, whereas a Poissonian PND is taken to a 
super-Poissonian PND, all moving towards the fixed thermal PND if re < 1. 



IV. Entanglement breaking property 

It is known that the phase conjugating channel is entanglement breaking [lH, It is also known that every 
entanglement breaking channel has a description in terms of rank one Kraus operators (35j . We demonstrate these 
aspects using our Kraus operators {Te(n)}. 

The Kraus operators Ti(k) presented in (|3.11j) are not of unit rank; indeed, rank Ti(k) = £ + 1, £ = 0, 1, 2, • ■ ■ . 
We noted immediately following (|3.11[) that T((k) are trace-orthogonal. In the generic case, trace-orthogonality 
requirement would render the Kraus operators unique, but this is not true with the present situation. The reason is 
that all these trace-orthogonal Tg (re)'s have the same Frobenius norm: tr to(/t)T^(«;)n = (1 + k 2 ) -1 , independent of t. 
Thus the set {T r } defined through T r (re) = {/^ r T^(re), for any unitary matrix (Ui r ) will be a set of trace-orthogonal 
Kraus operators describing the same channel as the original trace-orthogonal set {T^(re)}. 

More generally, and independent of trace-orthogonality, the map p — >• p = Eo^aWPVW describes the same 
channel as p — > p = J2gTg(K)pTj (re) if the matrix U connecting the sets {T^(re)} and {T a (n)} is an isometry [3(1 [H| : 

T a( K ) = E UtaTein), E u (aU* a = 5 ir 

a a 

J2 T ? ( k ) p T l ( k ) = E T « ( k ) p T ' 1 (») ■ ( 3 - 4 °) 

If the index set a is continuous, as in the case below, then is to be understood, of course, as an integral. Now, 
the matrix elements between coherent states |a) and Fock states \k) define such an isometry 

U ia ee (£\a) = cxp[-|a| 2 /2]^=. (3.41) 
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The resulting new Kraus operators T a {n) are 

2 °° I 



y/T\ 

1 e 

e ~^~ E (VT+^r 1 £ ^(vT + ^r n (VT+^=2)-^-")|^ - n)(n| 



=0 

of 



00 



£=0 n=0 

L=|a/i/i + K- 2 )K/Vl + k 2 |, Va e C. (3.42) 



\/T+ K 

It is manifest that rank T Q (ft) = 1 for all q€C, the complex plane, showing that the phase conjugation channel is 
indeed entanglement breaking. However {T q (k)} are not trace-orthogonal even though {Ti(n)} from which the former 
are constructed were trace-orthogonal. This is due to the fact that the isometry U defined in (|3.41j) is not an unitary, 
which in turn is a consequence of the overcompleteness of the coherent states. 

This brings us to another aspect of 2?(k). In terms of these new Kraus operators the phase conjugation channel 
T>(k) reads 



p =7T- 1 j d 2 ct T a (ft) p T a t (k) 

= 7r- 1 (l + ft 2 )- 1 / <PaQ((y/l + K 2 )-V)|a/Vl + K- 2 )(a/y/l + ft- 2 |. (3.43) 



Thus the diagonal weight function of the output state of the channel is the Q-distribution of the input state p : 



fout 



= k 2 Q m (n a*). We may combine this result with the earlier one on rank one Kraus operators to state 



Theorem 5 The diagonal weight of the output of the quantum-limited phase conjugation channel is essentially the 
Q -distribution of the input state. The channel T>(n) is not only classicality breaking, but also entanglement breaking. 

The diagonal weight of the output state at a is the Q-distribution of the input state evaluated at K _1 a*. 
Since Q(a) > for all a and for any p, the channel is nonclassicality breaking. The intimate relationship 
between this result and the earlier one on nonclassicality breaking may be noted. While the former followed di- 
rectly from the behaviour of the characteristic function, the present one required consideration of the Kraus operators. 



IV. BEAMSPLITTER/ATTENUATOR CHANNEL &(/«), < k < 1 

The two-mode unitary operator corresponding to the beamsplitter channel induces the following symplectic trans- 
formation on the quadrature operators of the bipartite phase space [l7j | : 



/cos6> -sinS 

cos 6 — sin ( 

sin 6 cosfl 

\ sin# cos( 



(4.1) 



Note that 5* is a direct sum of identical two-dimensional rotations: as in the case of T>(k), the position and momentum 
operators are not mixed by this transformation. The position variables transform as 



<l2 



M 



cos a sin V 
- sin 9 cos 9 



(4.2) 



and, consequently, the momentum variables as 



(Af 



-1\T I Pi 
P2 



M 



(4.3) 
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It is evident from S that the parameter k in Ci(k) is related to 9 through cos 9 = k, sin 6* = \/l — k 2 . The function 
F(z\, z 2 , rji, 772) of (|2.10p for the present case is given by 



F(z 1 ,Z 2 ,T] 1 ,r] 2 ) =exp 772 (\A - K 2 Zl + KZ 2 ) +77l(KZl - \f\ - k 2 z 2 ) 



(4.4) 



As in the previous case of T>(n), the differentiation on F(z-\ , 22, 771, 772) can be performed in a straight forward manner 
to obtain the matrix elements of the unitary operator |53j , leading to 



CI 



1 n 2 V n i \n2\m1\m2 ! 



E"i EZo ni a" 2 Q (-i)™ 2 -^' K^-^-j (Vi - « 2 ) 



r-\-ri2—j 



-\-7i2—r—j ■ 



(4.5) 



Now, to obtain the Kraus operators from these matrix elements we set, as in the case of T>(n), n 2 = and 777,2 = 
Setting 772 = =>■ j = 0, and we have 



= XI \/ m+£ Q ( - K 2 ) £ K m |m) (777 + ^|, £ = 0,1,2, 



(4.6) 



as the Kraus operators of the beamsplitter or quantum-limited attenuator channel. It is easy to see that the Kraus 
operators are real and pairwise trace-orthogonal, as in the case of T>(k). 



I. Properties of the Kraus operators 

We now explore the properties of the Kraus operators presented in Eq. (|4.6j) . connecting the action of the channel 
on the Fock basis to that on the characteristic function. We firstly exhibit the fact that the beamsplitter channel 
simply effects a scaling on the weight function of the diagonal representation. We show that vacuum is the only fixed 
point of the channel. It is further shown that in any set of Kraus operators describing the channel Ci(k), there will be 
not even one operator of unit rank, thus demonstrating that Ci(k) is not an entanglement breaking channel. Finally, 
the manifestation of the semigroup structure of the family of channels Ci(k), < k < 1 is brought out in the Kraus 
representation, as also an associated Zeno-likc effect. 

Recall that the beamsplitter channel induces the following transformation on the characteristic function fl7| : 

Xw(O^Xw'(0=Xw(nO exp[-(l - K 2 )|e| 2 /2] 

= *w(«0 cxp[ K 2 |C| 2 /2]exp[-|e| 2 /2]. (4.7) 

Thus the normal ordered characteristic function xat(£) transforms as 

X n(0 = Xw(0 exp(|C| 2 /2) -> Xn (0 - Xn(*0- (4-8) 

Since Xn(£.) and the diagonal weight (f>(a) form a Fourier transform pair, it is immediately seen that </>(«) gets simply 
scaled under the action of the Ci(ft) channel: 4>(a) — > </>' (a) = K~ 2 0(K -1 a) [54j |. 

It is instructive to bring out this fact from the perspective of the Kraus operators. Since every state p can be 
expressed through a diagonal 'weight' 4>(a) as [13] 

.-1 / J2 



p = n- L J d z acj)(a)\a){a\, (4.9) 

to exhibit the action of the channel on an arbitrary state it is sufficient to consider its action on a generic coherent 
state. We have 



\a){a\ ^J2B e (K)\a){a\Bl( K ) 



E E E ((1 "V |a|) (^r(^rS^HH ; (4.10) 
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where we used the fact that the operator 

oo 

\m)(n\ -^J2 B ^)\m)(n\Bl{K) 



1=0 

min{m,n} 



= J2 ^"^"^(l-K 2 ) K m+n - 2e \m-£)(n-e\. (4.11) 

1=0 

Carrying out the summations in Eq. (|4.10[) , one finds [5r| 

oo 

Y^B t {K)\a){a\B\{K) = \kol)(koi\. (4.12) 

1=0 

With this the action of the channel Ci(k) reads 

p — > p = 7T _1 / d 2 a 4>{a)\Ka) (kq| 



= tt -1 *;- 2 / d 2 a(j){K- 1 a)\a){a\, (4.13) 

which means 

Ci(k) : 0(a) -> /t- 2 ^ (k -1 ^ . (4.14) 
We have thus proved in the Kraus representation 

Theorem 6 The scaling <fi p (u) —¥ 4> p (ot) — K~ 2 (f> p (K~ 1 a), < K < 1, is a completely positive map whose Kraus 



decomposition is given by {B((k)} of 
As an immediate consequence we have 

Corollary 1 The beamsplitter channel cannot generate or destroy nonclassicality. 

Proof: By definition a state is classical if and only if its diagonal weight function (f>(a) is pointwise nonncgative 
everywhere in the complex plane [47| . Since a pointwise positive function goes to a pointwise positive function under 
the above scaling transformation, it follows that a classical state (and a classical state alone) is taken to a classical 
state under the action of the (quantum-limited) attenuator channel. 

II. Fixed Points 

We now examine, using cumulants, if there are any fixed points for this channel. In view of Eq. (|4.7[) . the cumulant 
generating function transforms as follows under the action of this channel : 

Ci(k) : r(0^r'(£)=log[ X vHO]-i(l-K 2 )|£| 2 (4-15) 

As in the previous case of T>(k), the cumulants of order > 2 of XwiO are 

/ d \ mi ( d 



7^ = («r i+ro2 ^j [o^) 

= (Kr i+ " i2 7 ™, (4-16) 

where 7 mim , and 7 mim2 are the cumulants of respectively the output and the input states. Thus, with the exception 
of the trivial case k = 1, action of the channel C\{n) attenuates the higher order cumulants. Therefore any state that 
is preserved is necessarily Gaussian. Since the additional term i(l — k 2 )|£| 2 in Eq. (|4.15j) does not involve a cross 
term between the real and imaginary parts of £, in looking for fixed points it is sufficient to consider the action only 
on thermal states. Given an input thermal state pth as in (|3.30[) the output, in view of Eq. (|4.1ip . is 

Pout = , { )r X \, Ef 1 ( f X 2 , V|n)(n|, (4.17) 
1 - (1 - n 2 )x ^ V 1 - (! - K W 
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where we have used the identity in Eq. (|3.32[) . Since p ut is a thermal state, comparing Eqs. f|3 . 30[) and (|4. 1 7[) we 
have this transformation law for the thermal parameter : 

2 

ft X 

Ci(re) : x -> — or a -> ft 2 a + (1 - ft 2 ). (4.18) 

1 — (1 — K Z )X 

This means that the output thermal state is always strictly 'cooler' than the input thermal state. Thus any thermal 
state is driven towards the ground state under repeated use of the channel. Indeed, an arbitrary state is driven 
towards the ground state. We thus have 

Theorem 7 For the action of the beamsplitter channel C\{k), the ground state is the only fixed point. 

As an illustration of the action of Ci(k), consider the case of a Fock state as the input. By Eq. f|4. 1 1[) we find 

n 

J2Be(n)\n)(n\Bl(n) = n Ct (1 - k 2 )V"~ 2 'I" ~ ( n ~ l \ 
t i=o 

n 

= ^"<^(1-K 2 ) n -V^)(£|. (4.19) 
e=o 

That is, the beamsplitter channel takes a Fock state to a convex sum of all Fock states with photon number less 

than or equal to n. It is thus clear that any input state which is diagonal in the Fock basis is taken to a Fock diagonal 
state at the output. We also note that for an arbitrary input p, a Fock diagonal entry \m)(m\ at the output gets 
contribution only from the Fock diagonal entries |n)(n| of the input with n > m. Putting these facts together, it is 
easy to see that no state which is in the support of a finite number of Fock states is preserved under the beamsplitter 
channel, for the strength of the highest Fock state strictly decreases. 

As a second example, consider the phase averaged coherent state (|3.38p as the input. By Eq. (|4.12[) . we have a 
phase averaged coherent state at the output : a Poissonian PND at the input is taken to another Poissonian PND at 
the output. Unlike the phase conjugation channel, both a thermal as well as a Poissonian PND at the input are taken 
at the output to respectively a thermal and Poissonian PND of strictly decreasing mean photon number. The fixed 
point should therefore be a Poissonian or thermal state of vanishing mean photon number; this is the vacuum state. 



III. The issue of Entanglement breaking 

It is known that the beamsplitter channel is not entanglement breaking [36|. It should thus be possible, as it is 
obligatory, to demonstrate that this channel cannot be represented using a set of rank one Kraus operators. We begin 
by noting that in the limiting case k = 0, all our Kraus operators Bi(0) are of rank one. Indeed, (B((0)) mn = 5 m o5 n i- 
This singular limit corresponds to the quantum-limited A\ channel which is known to be entanglement breaking. We 
consider therefore the nontrivial case k ^ 0. It is manifestly clear that rank B^(k) = oo for all £ (for k ^ 0). If we 
represent this channel using another set of Kraus operators {-B r (ft)}, then these new operators should necessarily be 
in the support of the set of operators {Be(n)}. Thus a necessary condition that one is able to represent the channel 
{Bi («)} using rank one Kraus operators is that there be (sufficient number of) rank one operators in the support of 
{Bi(n)}. It turns out that there is not even one rank one operator in this support. Indeed, a much stronger result is 
true. 

Theorem 8 : There exists no finite rank operator in the support of the set {B((n)}, k 7^ 0. 

Proof follows immediately from the structure of the £^(k)'s : Bq(k) is diagonal, and the ran" 1 entry of Bi(k) is 
nonzero iff n = m + I. Any matrix in the linear span of {Be(n)} is of the form M = ^iCgBifa), and is upper 
diagonal. Let N be the smallest I for which the c-number coefficient a ^ 0. Let M be the matrix obtained from the 
upper-diagonal M by deleting the first columns. Clearly, rank M = rank M. Further, the diagonal entries of the 
upper triangular M are all nonzero, being the nonzero entries of Bn{h). Now, the rank of an upper triangular matrix 
is not less than that of its diagonal part. Thus, rank M is not less than rank Bn(k) = 00, thus completing the proof. 
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IV. Semigroup property 

It is clear from (|4.7[) that successive actions of two beamsplitter channels with parameter values kx,k 2 is a single 
beamsplitter channel whose parameter k equals the product kik 2 of the individual channel parameters : 

Ci(ki) : xw{0 x'w(0=Xw(kiO exp [-(1 - k?)|£| 2 /2], 
Ci(«2) : xV«) -> x'V(0 = xV(KaO exp [-(1 - K |)|£| 2 /2] 

= X^(«x«20 exp[-(l-«i«2)l^| a /2]. (4-20) 

It is instructive to see how this semigroup property emerges in the Kraus representation. Let {Be 1 (kx)} and \Bg 2 (k 2 )} 
be the Kraus operators of the two channels. The product of two Kraus operators Bg 1 (rti), Bi 2 (fc 2 ), one from each set, 
is 

x j m + l i+ l *C tl+ t 2 {^2)™ K i 2 }\m){m + l 1 + £ 2 \. (4.21) 
Thus the action of the product channel on the input operator \r)(r + S\ is 
2 B* (ki)S^ 2 (K2)|r) (r + 5|B £2 ( K2 ) f B* («i) f 



x (V^ff (\A^)' 2 \l n+tl+ta c*i+t2 («i«2) n 4 

x | m )(m + ^+^ 2 |r)(r + (y|n + ^+^2>(n|. (4.22) 



Denoting ^1+^2 = ^, the expression on the RHS becomes 

r £ 00 



RHS =EE E '^^(l - *?f(i - Ki)^- <iJ (/ tl « 2 r+ n 

£=0 £i=0 m,n=0 



x y/t+^Ct^Ct S r>m+e S r+s , n +e \m)(n\. (4.23) 

The sum over l\ is the binomial expansion of [(1 — fCx) K 2 + (1 — = (1 ~ an0 ^ m addition, we have the 

constraints m + £ = r and n + i = r + 6. With this the expression (|4.23[) reduces to 



RHS ^^{l-K\Kl) l ^C^C' l {K l K 2 ) 2r - 2i+s \r-£){r~£ + 8\. (4.24) 



£=0 



Comparing Eqs. (|4.11[) and (|4.24[) we find that the expression in (|4.24[) is precisely the action of a quantum-limited 
attenuator channel with parameter k\k 2 : 

00 00 

B^ 1 ( Kl )B^( K2 )|r)(r + ^|B]j« 2 )B] i ( Kl ) = ^B^( KlK2 )|r)(r + <5|B]( Kl / t2 ). (4.25) 
ti,Zv=0 1=0 

An identical result can be similarly obtained for the behaviour of \r + 5)(r\, and thus we have proved the semigroup 
property 

Cx(kx) oCi(k 2 ) = Ci(«i«2). (4.26) 
Remark on interrupted evolution and Zeno-like effect : 

As seen from (|4.ip the parameter k specifying the channel Ci(k) equals cos 8, where 9 is a measure of the two-mode 
rotation effected by the 'beamsplitter' coupling the system mode to an ancilla mode assumed to be in the vacuum 
state initially. The associated two-mode unitary operator U{9) = cxp[— 9(a i b — b^a)} may be viewed as effecting 
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FIG. 2: Showing the variation of the attenuation factor k as a function of the evolution parameter 9, with the evolution 
interrupted periodically once every N~ n/2 in 9. The broken curve represents the uninterrupted evolution for < 9 < iv/2, 
total attenuation (k = 0) being achieved at 9 — n/2 in this case. The other curves correspond to N = 2, 3, 5, and 10. The Zeno 
tendency of k, for large N, to become linear in 9, with slope ~ iV -1 may be noted. 



evolution for 'duration' 6 under the Hamiltonian — i(afb — b^a). It is clear that attenuation increases monotonically 
as 9 varies from to tt/2, with total attenuation achieved at 9 = ir/2. 

Two-mode evolution for duration N~ 1 ir/2 followed by tracing away of the ancilla mode results in the channel 
Ci(kn,i) where Kn,i = cos(N~ 1 w/2). Now suppose that we have interrupted evolution in the sense that this process 
leading to Ci(ftjv,i) is repeated £ times. By the semigroup property the net result is a quantum-limited attenuator 
C\{hn,i) where knj = [i^nsY = (cos(A r ~ 1 7r/2))^. The behaviour of the attenuation factor knj is depicted in Fig. [2] 
That the effect of interruption is to slow down attenuation is transparent. For large N we have km,i ~ 1 — jjy (-4) 
reminiscent of quantum Zeno effect [5o| . 



V. AMPLIFIER CHANNEL C 2 {k), k > 1 

The two-mode metaplectic unitary operator describing a single-mode quantum-limited amplifier channel corresponds 
to the following symplectic transformation on the mode operators [j~7j : 

f'coshi' sinhi/ \ 

g cosh v — sinh^j f5 1) 

' sinh^ coshi/ I \ • > 

— sinh v cosh v ) 

As in the earlier two cases of T>(k) and Ci(k), the position and momentum variables do not mix under the action of 
C2{k). The position variables transform as 

q A^( q i)=M ( qi ) = ( COsh " M (5 . 2) 

?2/ \q 2 J \1 2 J \-svahu coshzy J \q2j 
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and the momentum variables transform according to M 1 . Thus the parameter re in C2(re) is related to the two-mode 
squeeze parameter v through re = cosIiza The function F {z\, z 2 , ?7i, 772) m (|2. 10[) is readily computed to be 

F(zi, «2, 771,172) = «;~ 1 exp|/t~ 1 (?7i2;i + r) 2 z 2 ) + ( v 7 ! - «~ 2 )( 7 7i ? 72 - 2122)} • (5.3) 

As in the earlier cases of T>(n) and Ci(re), the differentiation on F(z\, z 2 , fji, t? 2 ) can be performed to obtain the matrix 
elements of the unitary operator corresponding to the symplectic S in (|5.1|) . We obtain, after some algebra patterned 
after the earlier two cases, 

t- 1 

C™ 1 ™ 2 = = m!m 2 ! 

Vni!ri2!mi!wi2! 

x g ^*^i C . (_l)r (^1 _ „-2)r+m 1 -j (fc-ljn^-r^^ ^^^^ _ (5.4) 
r=0 j'=0 

The Kraus operators are obtained from C™ 1 ^ 2 by setting ri2 = 0, and m 2 = ^. Setting n 2 = =>■ r = 0, and we have 
A e (n) = K - 1 J2 V m+e C e (Vl - k- 2 ) (K-'nm + ^ml, * = 0, 1, 2, ■ • • (5.5) 

ro=0 

as the Kraus operators of the quantum-limited amplifier channel C 2 (re), re > 1 [42j . 



I. Duality between the attenuator family Ci(-) and the amplifier family C2( ) 

The Kraus operators Ae(re), re > 1 of the amplifier channel C 2 (re) have an interesting dual relationship to the 
Kraus operators -Bf(re _1 ), re > 1 of the attenuator channel Ci(re~ 1 ). While X^fco Ae( K )^( K ) = 1 5 « > 1 an d 
Sfco^l( K ) = ■"■! K < 1' consistent with the trace-preserving property of C 2 (re) and Ci(re ), we have 

00 

^A,(re)4(re) = K - 2 1, 

«=o 
00 

5^B/(k')bJ(k') = (k')- 2 1. (5.6) 

£=0 

Thus the (trace-preserving) families Ci and C 2 are not unital. But they are 'almost unital', for the failure to be unital 

is by just a scalar factor. This shows that the family {reA^(re)^, re > 1} and the family {re B((n )t, re < 1} too 
describe trace-preserving CP maps, and we may ask what these 'new' channels stand for. 

The meaning of these channels may be easily seen by considering the adjoints Ag(Ky, re > 1 of the Kraus operators 
of the amplifier channel : 

A e (^ = re" 1 J2 V m+e C e (V 1 - re~ 2 ) K- m \m)(m + £\ 

m=0 

= re-^re" 1 ) (5.7) 

Thus {reAf(re)^}, re > 1 are the Kraus operators of the beamsplitter channel Ci(re ) with re = re -1 < 1. Similarly it 
can be seen that {re Bg(K re < 1 represents the amplifier channel C 2 (re) with re = (re ) _1 > 1. Thus we have 

Theorem 9 The amplifier family C2(re) and the attenuator family Ci(re _1 ), re > 1 are mutually dual: their Kraus 
operators are connected through the adjoint operation. 



II. Properties of the Kraus operators 



We now explore the properties of the Kraus operators of C 2 (re) exhibited in Eq. (|5.5p . relating the action of these 
operators in the Fock basis to the defining or expected transformation property of the quasiprobability in phase 
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space. Wc begin by establishing that the amplifier channel simply scales the Q function. Then we show, through 
consideration of the cumulants, that there is no fixed point for the amplifier channel. The structure of the Kraus 
operators is invoked to comment on the fact that the quantum-limited amplifier channel is not an entanglement 
breaking channel. Finally, the semigroup structure of the quantum-limited amplifier channels is brought out in the 
Kraus representation. 

Under the action of the amplifier channel C2{k) the Weyl-ordered characteristic function transforms as follows, and 
this may be identified with the very definition of the channel : 

Xw(Z) -> Xw(0 = exp [-(k 2 - l)|£| 2 /2]. (5.8) 

Given a Wcyl-ordcrcd characteristic function xw(£)> the corresponding antinormal ordered characteristic function 
corresponding to the Q-distribution is [45| 

XA(0=Xw(0 exp[-|e| 2 /2]- (5.9) 
Therefore the channel action Eq. (|5.8p . written in terms of Xa(£,), reads 

xa(0->xA(0 = xa(«0- ( 5 - 10 ) 

That is, Xa(0 simply scales under the action of the amplifier channel, a fact that should be profitably compared with 
the scaling behaviour (|4.8)1 for the attenuator channel. Since xa(£) and the Q- function form a Fourier transform 
pair, the action of the amplifier channel is fully described as a scaling transformation of the Q-function: Q(a) — > 



Q'(a) = k~ 2 Q(k~ 1 o), k > 1 |H 



It is instructive to see in some detail how our Kraus operators Ai{k) bring out this behaviour. Given a state 

oo oo 

P= \n)(n\p\m)(m\ = pnm\n){m\, (5.11) 

n,m— n,m— 

its corresponding Q function is (45| 

°° (a*) n (a) m 

Q p (a) = (a\p\a) = exp[-|a| 2 ] V V > V > p nm . (5.12) 

— n Vni ymi 

n,m— U 

To see the action of the linear map C2(k) on an arbitrary p, it is sufficient to exhibit its action on the operators \n) (m|, 
for all n, m > 0. We have 

oo 

\n)(m\ -*^2A e (K)\n)(m\Al(K) 

1=0 

/ \ — (n+m) oo _ 2\£ 

= k- 2 [K) J2 ~l V^+W-Vi^+ey-\n + i){m + t\. (5.13) 

Vn.m. e=0 

Thus, under the action of the channel C2(k), p goes to 

oo —(n+m) 00 /i — 2\i 

P = K - 2 £ P^^TTT E t\ \^+W-V^+W-\n + l)(m + £\. (5.14) 

n.m—O v £—0 

The Q function of the resultant or output state p is 



(a\p'\a) = K - 2 exp[-H 2 ] £ ^^(aTW" 1 )£ \a\ M 

n,m=0 v \£=0 J 



k 2 exp[-| K ^l 2 ] V ^ " ^ ^ p nm 
— n \n\ vm! 

n,m— 



= K- 2 Q(K- 1 a). (5.15) 

We thus conclude 
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Theorem 10 The scaling Q p (a) — > Q > (a) = k 2 Q p (k a), < k 1 < 1, is a completely positive map whose Kraus 
decomposition is given by {Ai(k)}. 

This result may be compared with Theorem 6 for the Ci(-) family of channels. 

The amplifier channel has the following property in respect of nonclassicality of the output states : 

Corollary 2 The amplifier channel cannot generate nonclassicality. 

Proof: By Eq. (|5.8|) . the normal ordered characteristic function transforms as follows 

C 2 (k) : Xn(0 -> Xn(0 = Xw{<) exp [-(k 2 - 2)|£| 2 /2]. (5.16) 
This may be rewritten in the suggestive form 

Xn(0 -> X N (0 = XN«) cxph( K 2 - 1)|£| 2 ]. (5.17) 

Fourier transforming, we see that the diagonal weight <f> (a) of the output state is the convolution of the (scaled) 
input diagonal weight with a Gaussian (corresponding to the last factor), and hence it is pointwisc nonncgative 
whenever the input diagonal weight <p(a) is pointwise nonnegative. 

Remark : We are not claiming that the amplifier channel cannot destroy nonclassicality [compare the structure of 
Corollary 2 with that of Corollary 1 following Theorem 6]. Indeed, it is easy to show that nonclassicality of every 
Gaussian state will be destroyed by any Ci{K) with n > s/2 fill - 1431 l57ll . It is also easy to show that there are states 
whose nonclassicality will survive C2(k) even for arbitrarily large n [4ll. l42l l57j. To see this, note first of all, that any 
state p whose Q-function Q(a) = (a\p\a) vanishes for some a is necessarily nonclassical. The assertion simply follows 
from the fact that under the scaling Q(a) — > K _2 Q(K _1 a) a zero ao of Q(a) goes to a zero at Kao- 



III. Fixed points 

By Eq. (|5 .8[) we have, under the action of the channel C2(k), the following behaviour for the moment generating 
function : 

r(0 -> r'(o = log [ X w(kO] - i|CI 2 (^ 2 - 1). (5.18) 

The cumulants of order > 2 of the output characteristic function are 

= {K) m ^ lmim2 , (5.19) 

where r y mi m 2 are the cumulants of the input state. Thus, for any non-Gaussian input state, the higher order cumulants 
grow monotonically with repeated use of the channel. Thus, leaving out the case k = 1 which corresponds to the 
identity channel, there is no non-Gaussian state that is preserved. To see if there is any fixed point among the 
Gaussian states, it is sufficient to consider only thermal states as input (|3.30l) . By Eq. (|5.13[) . the output state is 

oo 

Pout = (1 - x)k~ 2 J2 k~ 2i V - 1 + x) n \n){n\, (5.20) 

n=0 

where we used the binomial expansion to perform a sum. Comparing Eqs. (|3.30[) and (|5.20p . we have 

C2(ft) : x — > x = k~ 2 (k 2 — 1 + x) or ao — > a = n 2 ao + k 2 — 1. (5-21) 

It is clear that the output thermal parameter a is strictly greater than ao- Hence there is no thermal state that is 
left invariant under the action of the amplifier channel. Collecting these facts together, we conclude 

Theorem 11 There exists no state which is a fixed point of the guantum-limited amplifier channel. 
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As a simple illustration, consider the action of the channel on a Fock state. We have by Eq. (|5,13l) 

oo 

C 2 (/s) : \n)(n\^ K- 2 J2 n+e Ce(l-K-yK- 2n \n + e)(n + l\. (5.22) 
e=o 

Thus an input Fock state |n)(n| is taken to a convex combination of all Fock states with photon number greater than 
or equal to n. It may be noted that this behaviour is complementary to that of the beamsplitter channel where the 
output was a convex combination of all Fock states upto \n)(n\. And wc find, analogous to the beamsplitter case, that 
any state which is in the support of a finite number of Fock states is not preserved by the channel. 

As a second example, consider the phase averaged coherent state (|3.38[) as the input. By Eq. (|5.22l) the output is 

Pout = J2 j+iCi ~~~r~( 1 - K- 2 y K - 2 3\j+£)(j+£\. (5.23) 

It is clear that we cannot solve consistently for parameters k > 1 and A such that the output PND is also Poissonian. 
In view of Corollary (2] we may conclude that an input Poissonian PND generically results in a super- Poissonian 
PND. 

Remark on entanglement breaking : It is well known that the quantum-limited amplifier channel is not entan- 
glement breaking [361 ]. It may be pointed out in passing that this fact follows also from the structure of our Kraus 
operators {Ae(n)}. Since these operators coincide with the transpose of the beamsplitter channel Kraus operators 
{Be^n -1 )}, apart from a ^-independent multiplicative factor, there exists no finite rank operator in the support of the 
set of operators {Az{k)}. In particular, there are no rank one operators in the support of {Ai(k)}. Hence, C 2 (k) is 
not an entanglement breaking channel. 



IV. Semigroup property 

It follows from the very definition of the amplifier channel that the composition of two quantum-limited amplifier 
channels with parameters K\ and n 2 is also such an amplifier channel with parameter k\k 2 > 1 : 

C 2 (k 2 ) o C 2 ( Kl ) : xw{0 -> Xw(0 = Xw(ki*2 exp - l)|£| 2 /2]. (5.24) 

That is, 

C 2 (n 2 ) o C 2 {ki) = C 2 (kik 2 ) = C 2 (ki) o C 2 (k 2 ). (5.25) 

It will be instructive to examine how this fact emerges from the structure of the Kraus operators. Let the set {Ai x («i)} 
be the Kraus operators of the first amplifier and let {Ae 2 (K2)} be that of the second. Then the product of a pair of 
Kraus operators, one from each set, is 



x ( KlK2 )-" Kl - £ > + 4 +£ 2 )(n\. (5.26) 
Thus, under the successive action of these two amplifier channels the operator + 6\ goes to 

E A ei (Ki)AtM\j)(j + 5\A t . 2 {K 2 y A ei ( Ki y 
ti,t 3 

OO OO 

= (*i*>r a EEE ll+l2C ^ (i - -i 2 t (i - «?t («i- 2 )- (n+m) -r % 

i!,i 2 n=0 m=0 

x c tl+l2 m +^+ e - C ei+ta \n + h+ i 2 ) (n\j) (j + S\m) (m + £ x + £ 2 \ . (5.27) 

Denoting l\ + £2 = £, the right hand side of the above expression reduces to 

00 £ 00 00 

(^r 2 EEEE ^ (i - -r 2 ) £i 2 a - - 2 2 )r^ ] M~ {n+m) 

=0 ^i=0 n=Q m=0 



X y/"+lCt m+l C t 6 m ,j+8 S n ,j |n + £){fl + S + £\ (5.28) 
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FIG. 3: Showing the evolution of the amplification factor k as a function of the parameter [i, with the evolution interrupted 
every 2TV -1 in /i. The broken curve represents uninterrupted evolution for < /i < 2. The other curves correspond to 
TV — 2, 3, 5 and 10, the total duration adding to 2 for each TV. The Zeno tendency of k, for large TV, to become linear in fj, with 
slope ~ TV" 1 should be noted. 

The summation over the index l\ is a binomial expansion : 

£ *c tl (i - «r 2 ) £l («r 2 (i - ^ 2 )f- ei) = (i - «r W- 

Thus the expression in (|5.28[) reduces to 

oo 

( KlK2 )- 2 J2 (1 - «r 2 ^ 2 )' («iK2)- W+J+tf) 0' +/ C<J+'+ s C< I j + 00' + * + <*|- (5-29) 

Comparing Eqs. f|5 . 13[) and (|5.29[) . we see that the latter is the Kraus representation for a single quantum- limited 
amplifier channel. That is, 

Y,AiMAiM\3){3+S\A l2 {K2) ] At,{Ki)^ =Y J M^2)\j){j + S\Ai{K l K^. (5.30) 

l x ,li £ 

A similar behaviour holds for \j + 5}(j\ as well. And this is what we set out to demonstrate. 
Remark on interrupted evolution and Zeno-like effect : 

As seen from (|5.ip the parameter k of C%{k) equals cosh/i, where /z is a measure of the two-mode squeezing effected 
with the help of an ancilla initially in the vacuum state. The relevant two-mode unitary squeeze operator U{fx) = 
exp[— /x(a^ — ab)] can be viewed as evolution for a duration /j, under the Hamiltonian — i(a'b' — ab), the amplification 
factor k = cosh \x increasing monotonically with increasing (i. 

For convenience let us fix the total duration in /i to some value, say 2. In place of this single evolution, let us 
consider a sequence of TV interrupted evolutions, by tracing away the ancilla after every 2TV" 1 duration in /x, the total 
duration adding upto 2. The semigroup property of these interrupted evolutions after a total duration of 2 in ji will 
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be a quantum-limited amplifier C2(k), with k = [cosh (2/^)]^, which should be compared with k = cosh 2. The effect 
of interruption in slowing down amplification is clear. The behaviour of this interrupted amplification is shown in 
Fig. [3l After I such evolutions the amplification factor will be k = (cosh2/iV)^, which for large N has the behaviour 
(t«l|| (^r)i reminiscent of quantum Zeno effect (55j as in the case of quantum- limited attenuation in Section IV. 



VI. THE SINGULAR CASE Ai 



We now consider briefly A2, the last of the quantum limited Bosonic Gaussian channels. The two-mode mctaplcctic 
unitary operator representing A2 produces a symplectic transformation on the quadrature variables which does not 
mix the position variables with the momentum variables [I?} : 

9i\ ^ ( Qx \ - M I 



q-2j \q 2 J V? 2 

M =(;_\). (a.. 

Therefore, our general scheme applies to this case as well. Unlike in the earlier cases of T>(k), Ci(re), and C?,{k) : in the 
present case it turns out to be more convenient to evaluate the matrix elements of J7'° & ' in a mixed basis : 

C%& = (mi\WW\n 1 )\n i ). (6.2) 

Here \q) labels the position basis of the ancilla mode. With this mixed choice, the Kraus operators are labelled by a 
continuous index 'q\ and are given by 

V q = (q\uW\0) = C/o'l^iXml. (6-3) 

mi ,ni 

where 

= dqi{m 1 \qi){q,q 1 -q\ni,0). (6.4) 



Here we have used, as in the earlier cases, the action of the unitary operator in the position eigenstates of the two-mode 
system. Employing the position space wavefunctions of the Fock states, we have 



-3/4 

rmiq 



CCo' - " mi , H ni (g)e-^ / d qi H mi ( qi )e-^-^. (6.5) 
^/2" 1 + mi ni!m 1 ! J 

The above integral is easily evaluated [58|, and we have 

C Zo q = , "'I 7 ', , Q mi H ni (q)<xp[-tf/4] 
v / 2" 1 + TOl ni!mi! 

= (m a \q/V2)(q\ ni ), (6.6) 

where \q/y/2) is the coherent state \a) for a = q/y/2, and the purpose of the round bracket being to distinguish the 
same from the position eigenket \q/y/2). With this notation the Kraus operators are 

V q = \q/V2){q\. (6.7) 

That the trace-preserving condition on the Kraus operators is satisfied emerges from the fact that the position kets 
are complete: / dqV^V q = J dq\q)(q\ = 1 . 
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To connect these Kraus operators V q to the action of the channel in the phase space picture, we examine the 
behaviour of an arbitrary pure state under passage through the channel. We have 



A 2 : p = IVXVI ->• P = I dq \q/V2) (?|l&><1%) 

dq mq)\ 2 \q/V2) (q/V2\ 



dqdp \^(q)\ 2 8{p)\[q + ip]/V2) ([q + ip]/V2\. (6.8) 

The last expression is already in the 'diagonal' form in the coherent states basis, with \ip{q)\ 2 8{p) 1 a = (q + ip)/v2 
forming the diagonal weight function 4>(ot). It follows by convexity that for an arbitrary input state p the output of 
the channel is given by 

'<Pa<Ka)\a)(a\, <f>(a) = (q\p\q) 6(p). (6.9) 

It is seen that this transformation is the same as xw{C) ~^ Xw ( < - 1H 2 ' T3 - > £j cxp[— 1£| 2 /2], the expected behaviour of 

the characteristic function under passage through _4 2 HI] . 

The above results can be alternatively understood through the action of the channel in the Fock basis. Under 
passage through the channel, 

\n){m\ -> J dqV q \n){m\V^ 

= J dq\q/V2) (q\n)(m\q) (q/V2\ 

— 1/2 

dq ^ H n (q)H m (q) \q/y/2) (q/V2\, (6.10) 

V2"+ m n!m! 

for all n, m. The outcome for an arbitrary input state p follows by linearity, and we have 
Theorem 12 The channel A2 is both nonclassicality breaking and entanglement breaking. 

Proof: We note from Eq. (|6.7[) that the Kraus operators are already in rank one form, thereby showing that the 
channel is entanglement breaking. And from Eq. (|6.9p we see that the output of the channel, for every input state p, 
supports a diagonal representation with nonnegative weight (g|p|(7) S(p) > 0, for all a = (q + ip) /y/2, showing that 
the output is classical for all input states. 

Remark on fixed points : It can be seen directly from the action of the channel on the characteristic function that 
the momentum variable of the output is set to zero, and then multiplied by a Gaussian in that variable. Thus any 
cumulant of order > 2 in this variable is set to zero by the channel action. Moreover it is easily seen that the second 
moments are not preserved, thus showing that there is no state that is invariant under the action of this channel. 



VII. SINGLE QUADRATURE CLASSICAL NOISE CHANNEL Bi(a),a>0 

The channel B\(a) 1 whose action is to simply inject Gaussian noise of magnitude a into one quadrature of the 
oscillator, is neither quantum-limited nor extremal. It can be realized in the form 

Bi(a) : p^p = -L / dqexp[-q 2 /a]D(q/V2)pD(q/V2^, (7.1) 
V7ra J 

where -D(a)'s are the unitary displacement operators. B\(a) is thus a case of the so-called random unitary channels 
[60j . a convex sum of unitary channels. The continuum 

Z q = (nay 1/4 cxp[- 2 /2a] D(q/V2) (7.2) 

are the Kraus operators of this realization. The quantum-limited end of Si (a) is obviously the identity channel, 
corresponding to a — > [lim a _5.o V™ exp[— q 2 /a] = 5(q), and Z q=0 = identity]. One may assume a = 1 without loss 
of generality. 
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The reason we present a brief treatment of this channel here is just to demonstrate that this case too subjects 
itself to our general scheme presented in Section II. A further reason why we treat this noisy channel here is this. In 
Section IX we shall treat every noisy channel as the composite of two quantum-limited channels, the case of Bi(a) 
constituting the only exception wherein this cannot be done. 

The two-mode metaplectic unitary operator representing B\ produces a symplectic transformation on the quadrature 
variables which, as in the earlier cases of T>(n), Ci(k), C2(k), and A2, does not mix the position variables with the 
momentum variables 1171 : 



M 



Qi 

1 -1 

1 



M 



And Pi,P2 transform according to (M 1 ) T . 

As in the immediate previous case A2, the matrix elements of [/'° & ' are 

{m x \{q\U^\n x )\n 2 ), 



(7.3) 



(7.4) 



where |g)'s are the position eigenvectors. In view of this the Kraus operators are labelled by a continuous index 'g' 
and are given by 



{<1 



(7.5) 



where 



C£tf = / dq 1 {m 1 \q 1 )(q 1 \{q\U^\n 1 )\0) 
dq 1 (m 1 \q 1 )(q 1 - q,q\m,0). 



Here we made the two- mode metaplectic unitary operator act on the position basis. 
To evaluate the Kraus operator, it is sufficient to evaluate the matrix elements 



C 



mq 
miO 



r-3/4 



V2 ni + mi f« 1 !m 1 ! 



dq 1 H„ 1 (q 1 - q)H mi (q 1 )e 2 e 2 



The above integral can be readily performed [611 ] . and we obtain 

-1/4, 



e 2 



e 4 



mi! 



nil VV2 

STr-^e-^^ilDCg/v^JIm) = Z q . 
We have thus recovered (|7.2I) . but staying entirely within our general scheme 



j m— mi („2 

mi 



(9 2 /2) 



(7.6) 



(7.7) 



(7.8) 



VIII. EXTREMAL GAUSSIAN CHANNELS 



It follows from the very identification of channels with trace-preserving CP maps that channels form a convex set, 
and since a convex set is fully characterized by its extremal elements, it is of interest to know if there are any extremal 
elements among the quantum-limited Bosonic Gaussian channels. The present Section is devoted to this issue. 

A convenient characterization of extremality in terms of Kraus representation is due to Choi [30| : A trace-preserving 
CP map f2 is extremal if and only if 57 supports a Kraus representation p — > p = J2 n W n p W\ , J2 n W n = 1, such 
that the set of operators {W^ W n }, with m and n independently running over their range, are linearly independent. 
While Choi's original result was formulated in the finite dimensional case, it is to be expected that this result will 
generalize to the infinite-dimensional case with suitable technical regularization. That this is the case is indicated, 
for instance, by Theorem 2.4 of the work of Tsui (62j . Since we are dealing with concrete physical situations, we may 
assume in what follows that such technical requirements are indeed satisfied. 

With these preliminary remarks in place, we are now ready to test the quantum-limited Gaussian channels for 
extremality. We begin with the transpose or phase conjugation channel. 
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I. The phase conjugation channel T>(k) 

The Kraus operators {T m (K)} for V{k) have been presented in p. lip . We may readily express the products 
{T^k) T n («;)}, < m,n < oo in the convenient form 

oo 

T} +s (K)T t (K) = Y,d(Sh\j + S)(j\, 

3=0 

OO 

r;( K )r £+5 ( K ) = X)d(<y)«li>0' + *l, 

3=0 

d(<*)« = [i + k 2 }-'- 5 ' 2 - 1 [(i + K- a rY " j for i < * 

= for j >i. (8.1) 

The two cases involving |^ + <5)(^| and + S\, S = 0, 1,2, • • ■ correspond respectively to m - > and m — n < 0. 
It is clear from DEHl that the set {T^(k) T„(k)}, < to, n < oo breaks into nonintersecting subsets labelled by 
m—n = 0, ±1, ±2 • • • . The index I = min{m, n} labels the TJj(k) T„(k)'s within a subset. That subsets corresponding 
to two different values of m — n are mutually orthogonal is manifest. Thus, we are left to examine only linear 
independence within each subset, and this is easily accomplished as follows. 

For a given value of to — n (i.e. fixed value of 5) arrange the c-numbcr coefficients d(S)ij into a matrix d(S) = 
((d(S)ij)). This matrix is lower diagonal, and its diagonals are nonzero [indeed, all the entries on and below the 
diagonal are non zero]. Thus d(6) is nonsingular. Therefore the linear independence of Tj (k)Tq(k), T^, s (k)Ti (k), 
T2 +(5 (k)T2(k), ■ • • follows from the linear independence of the mutually orthogonal |5)(0|, |1 + <5)(1|, |2 + <5}(2|, • ■ ■ and, 

similarly, the linear independence of Tg (k)Ts(k), T^(k)Ti + s{k), T\ (k)T2+5(k), ■ ■ ■ follows from the linear independence 
of|0)R |1)(1 + 5|, |2)(2 + <5|, .... 

This completes proof of the linear independence of the set {T^k) T„(k)}, < m,n < oo and hence proves 
cxtremality of T>(k) in the sense of Choi, for all k ^ 1. 

Remark on the Doubly Stochastic Case T>(1) : The channel T>(\) is exceptional within the phase conjugation 
family of quantum-limited channels in that it is both trace-preserving and unital, i.e. it is doubly stochastic. While all 
random unitary channels (convex sums of unitary channels) [63| are manifestly doubly stochastic, the converse is not 
true. Early (finite-dimensional) counter ex amp les to the converse can be found in Ref. (52l . |64| . It is this phenomenon 
that underlies a conjecture of Winter et al [38(. The channel T>{1) is a counter example from the infinite-dimensional 
Gaussian domain, and the only doubly stochastic case among quantum-limited Gaussian channels. That T>(1) is not 
random unitary readily follows from its cxtremality which we have established above. Not every non-extremal doubly 
stochastic map is random unitary, but an extremal doubly stochastic map can have not even one unitary operator in 
the support of its Kraus operators, for if it had it will be a convex sum of that unitary channel and another doubly 
stochastic map. The noisy classical noise channels of families B\ and B2 are obviously doubly stochastic and obviously 
random unitary. 



II. The attenuation/beamsplitter channel Ci(/t) 

The beamsplitter channel C\{k) is described by the Kraus operators {B m (K)} given in (14.61) . As in the earlier case 
of T>(k) we can compute the products {B^k) B„(/t)} in the form 

00 

OO 

b\{k)B 1+s {k) = + 

= for j < t. (8.2) 
Again the two cases \ j + 5}(j\ and \ + S\ correspond respectively to m — n > and to — n < 0. 
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The situation is similar to the earlier case of V(k). The set {B^^B^k)}, < m, n < oo fib-rates into nonintcr- 
secting subsets, and these are labelled by m — n = 0, ±1, • • • . The index £ = min{m,n} acts as the label within a 
given subset. The subsets being mutually orthogonal, it only remains to examine linear independence within each 
subset. 

As in the case of T>(k) we arrange the c-number coefficients f(8)ij into matrices f(8) = ((f(6)ej)), one matrix 
for each subset 8. These matrices are upper diagonal. None of the diagonal elements vanishes, and so the matrices 
are nonsingular, proving the linear independence of the sets {B\ +& (k)Bz(k)}i and {b\(k)Bi + s(k)\ t and hence the 
extremality of Ci(k), for all K > 0. 

III. The amplifier channel &(&) 

The amplifier channel C2(k) described by the Kraus operators {A m (n)} presented in (|5.5p turns out to be a little 
more subtle in respect of our present purpose. There is considerable similarity with the two earlier cases, though. As 
in the case of T>(k) and Ci(k), let us begin by expressing the products {A^k) A n (n)} in the form 



4 +s (K)A e (K) = jr,h(S) tj \j)(j + 6\, 



3=0 



A\(k)A, +s (k) = X>(%|j + 5){j\, 



3=0 



h(S) £j = K- 2 ^+j+ s C l+5 *+j+sc j+5 (^/i-T^ft+s 



(8.3) 



The two cases correspond respectively to m — n > and m — n < 0. The set of operators {A^^A^k)} manifestly 
separate into nonintcrsecting subsets, determined by m — n; the index £ acts as a label within each subset; the subsets 
are mutually orthogonal; and it only remains to determine the linear independence within each subset. Up to this 
point the situation is similar to the earlier two cases. 

As in the earlier two cases, let us arrange the c-number coefficients h(8)tj into matrices h(5), one matrix for each £. 
And now arises the distinction : whereas the invertibility of d(5) and f(5) was manifest, being lower or upper diagonal 
matrices, this is not so in respect of the present case. We are therefore led to demonstrate the nonsingularity of h(S) 
in a somewhat different manner. 



The multiplicative scalar (k 



) can be dropped from h(S) for our present purpose. We note further 



that h(5) can be simplified by left and right multiplication by diagonal matrices L, M : 



L , 



h(S) = LhR 

:S rS: Rrs — 



-2r 



y/r\(r + S)\ 



(8.4) 



Since L, R are nonsingular, invertibility of h(5) is equivalent to that of h(S). The new matrix h is symmetric and 
has a simple structure : h(6) rs = (r + s + 5)1. To demonstrate the nonsingularity of h(S), we begin by writing out its 
entries in detail : 



/ SI 
(6+1)1 
(6 + 2)1 



h(6) 



(6 + 1)1 
(6 + 2)1 
(6 + 3)\ 



(6 + 2)1 
(6 + 3)1 
((5 + 4)! 



(8 + l)\ 



(5 + k)\ 



(8 + k + £)\ 



\ 



\ 



(8.5) 



Dividing the first row of h(8) by 81, the second row by (8 + 1)!, the third row by (i5 + 2)!, and so on, we obtain 

(l (5 + 1) (8 + l)(6 + 2) (5+l)(5 + 2)(5 + Z) \ 

1 (5 + 2) (5 + 2)(5 + i) (6 + 2)(5 + 3)(5 + A) 

1 (6 + 3) (6 + 3)(5 + A) (8 + 3)(6 + A)(6 + 5) (8-6) 
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Now performing the row transformations R± — > Ri — i?o, R2 



(\ {5 + 1) (S+l)(5 + 2) (S+l)(8 + 2)(5 + 3) 
1 2(6 + 2) 3(<5 + 2)(<5 + 3) 

1 2(5 + 3) 3(8 + 3)(S + 4) 

1 2(5 + 4) 3(<5 + 4)(5 + 5) 



R2 — R% , and so forth we obtain 



(8.7) 



/ 



We can repeatedly do this kind of row transformations, starting with row R2 in the next iteration and row R3 in the 
subsequent iteration and so on. The matrix h(S) gets finally transformed to the form 



Z 1 









1 










2! ••• 








3! 











I : 







4! 



(8.8) 



which is upper triangular with non- vanishing diagonal entries, and thus has 'full rank'. We may thus conclude that 
the coefficient matrix h(S) is invertible, showing that the set {^4j +( j^4f }f [as well as {A^f+alf] is linearly independent, 
for each 5. 

This completes proof of linear independence of the set {A^k) A n (K,)}, < m, n < 00. Extremality of the amplifier 
channel C2(re) is thus established,, for all k > 1. 



IV. The Singular case A2 
Computation of the product v\ V q is nearly trivial in this case. We have from (|6.7[) : 

VW q = (q'/V2\q/V2)\q')(q\. (8.9) 

The inner product (q / \[2\qj V2) = exp[— q 2 /A — q /4 + qq /2] is nonzero for arbitrary q, q . Since \q )(q\ form a 
linearly independent set of operators, we conclude that the channel Ai is extremal. 

We have considered in this Section the quantum-limited channels 2?(ft), C\(k), C2(k), and Ai in that order. Since 
these are the only quantum-limited Bosonic Gaussian channels, the main conclusion of our study can be stated in the 
following concise form : 

Theorem 13 All quantum-limited Bosonic Gaussian channels are extremal. 

Remark on Ai : We have not treated separately the quantum- limited channel Ai, for this may be viewed as a 
particular case of Ci(re), corresponding to cos(#) = in (|4.1[) or, equivalently, k = 0. Thus, the assertion above covers 
Ai as well. 



IX. NOISY CHANNELS AS COMPOSITES OF QUANTUM-LIMITED CHANNELS 

Our considerations so far have been in respect of quantum-limited channels. We turn our attention now to the case 
of noisy channels. It turns out that every noisy channel, except B\(a) which corresponds to injection of classical noise 
in just one quadrature, can be realised (in a non-unique way) as composition of two quantum limited channels, so 
that the Kraus operators are products of those of the constituent quantum-limited channels. 

We have noted in Sections IV and V that composition of two quantum limited attenuator (or amplifier) channels is 
again a quantum-limited attenuator (or amplifier) channel. This special semigroup property however does not obtain 
under composition for other quantum-limited channels. In general, composition of two quantum-limited channels 
results in a channel with additional classical noise. For this reason we restore in this Section the original notation 
T>(k; 0), C\(k] 0), and Ci(k\ 0) in order to make room for this additional classical noise. 



30 



X,Y -)■ 
I 


23(«i;0) 


Ci(«i;0) 


C 2 (ki;0) 


-4a (0) 


V(k 2 ;0) 


Ci(«2«i; 2k 2 (1 + 
for < 1. 
C2(«a«i;2(l + k|)), 
for > 1. 


2?(K2/ti; 2re|(l — k| ) 


2?(/t 2 Ki; 0) 


A(2«l) 


Ci(k 2 ;0) 


2?(«2Ki; 0) 


Ci(k 2 «;i; 0) 


Ci(k2Ki;2«|(ki — 1)), 
for K2K1 < 1. 
C 2 (k 2 ki; 2(1 — «!)), 
for K2K1 > 1 


A(0) 


C 2 (k 2 ;0) 


2?(k 2 «i;2(^ - 1)) 


Ci(k 2 Ki;2(/s| - 1)), 
for K2K1 < 1- 
C 2 (k 2 ki; 2/tl(l — «?)), 
for K2K1 > 1. 


C2(«2«i;0) 


,4 2 (2(«|-1)) 


A(0) 


A(V K i+ 2 - 1) 


A(V2-«?-l) 


Aa(«i - 1) 





TABLE I: Showing the composition X o Y of quantum-limited channels X, V assumed to be in their respective canonical forms 
simultaneously. The composition results, in several cases, in noisy channels thereby enabling description of noisy Gaussian 
channels, including the classical noise channel 62(a), in terms of discrete sets of linearly independent Kraus operators. 



I. The composite C 2 (k 2 ;0) oCi(ki;0), k% > 1, Ki < 1 



It is clear from the very definition of these channels through their action on the characteristic function that the 
composite £2(^2; 0) o Ci(ki; 0) is a noisy amplifier or attenuator depending on the numerical value of ^2^1 : it equals 
C\(k2Ki; 2(«2 — 1)) for K2K1 < 1, and C^ij^i^i^ 2k|(1 — K i)) f° r K 2^i > 1, as may be readily read off from TablcUl 

The Kraus operators for the composite is given by the set {A m (K-i) B n {ni)} with m,n running independently over 
the range < m,n < 00. That these Kraus operators are linearly independent may be seen as follows. Since 
A m (K-2) is the same as B m (K 2 ~ 1 )^ except for an m-independent multiplicative constant, linear independence of the 
set {A m (K 2 )B n (Ki)} is the same as linear independence of {B^k^ 1 )^ B n (K,{)} . It will be recalled that in proving 
linear independence of the set {B m (K,y B n (n)} in Section V in the context of extremality of C2(k), we used only the 
structure of the -B^(-)'s in respect of the expansion coefficients in the basis {|m)(n|} being zero or nonzero, and not 
the actual numerical values of the nonzero coefficients. Hence the same argument should be expected to apply to 
proof of linear independence of {B^k^ 1 )^ B n (ni)} ~ {A m (K2)B n (Ki)} as well. That this is indeed the case is seen 
by computing the products A m (K 2 )B n (Ki) in the form 

00 

A £+s (K 2 )B t {Ki)=J2gi(S)i j \j + S)(j\, 
3=0 
00 

Ae(K2)B e+ s(Ki) = ^2gi(8)ej\j)(j + 6\, 
j=n 

g 1 (6) ij = n^ 1 V> +5 C e+5 W e (Jl - nj) (K2 1 K i y~ t , for j > I, 



= 0, for j < £; 



= 0, for j < I. (9.1) 

Comparing this with (|8.2[) . it is now clear that A m (K2)B n (KiYs are linearly independent for exactly the same reason 
by which B m (K)^ _B n («;)'s were linearly independent rendering the quantum-limited attenuator Ci(k;0) extremal. 
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II. The composite Ci(«2; 0) o £2(^1; 0), K2 < 1, fti > 1 

Again the composite Ci(k2]0) o C 2 (ki;0) is a noisy amplifier or attenuator depending on the numerical value of 
K2 K i, an d the details may be read off from Table |TJ The Kraus operators for the composite Ci(«2;0) oC 2 (ki;0) are 
given by {B m (K 2 )^4 n ( K i)} > < m,n < 00. To establish their linear independence we compute B m (n2)A n (Ki) in the 
form 

00 

B e+s { K2 )A e { Kl ) =J292(8)ej\j)(j +S\, 
3=0 

OO 

B/(/s a )^+«(«i) = £fla(*)«|j + <5)0"li 

<? 2 (% = «r x V^+^+^'+^cw '^ ^i-^r 2 )' ^ (V 1 - 4 

h(Sh = ^i 1 #V^ K i J (Jl - «r 2 ) + (V 1 - «l) ■ (9.2) 

Comparing with (|8.3[) we see that the Kraus operators B m (K2)A n (Ki) of C\[k2\ 0) o £2(^1; 0) are linearly independent 
for the same reason by which the quantum-limited amplifier channel was extremal. 




III. The composite T>(k2) oD(ki), «2,ki > 



Similar to the earlier two cases, the composite T>(k,2]0) o2?(k 1 ;0) is a noisy amplifier or attenuator depending on 
the numerical value of K2K1, and the details can be read off from Table HI It may be noted, again from Table H that 
this case tends to be more noisy that the earlier two cases. 

The Kraus operators for this composite are given by {T rn (K2)T n (ni)} , < m,n < oo. To exhibit the linear 
independence of these Kraus operators we compute the products T m (K 2 )T n (Ki) in the form 

oo 

Ii +4 (« 2 )r<(Ki) = ^2g 3 (5) e Aj + S)(j\, 

3=0 

oo 

2i(« 2 )Ii +a («i) = X>(fy 3 -|j + S)(j\, 



3=0 



93(S)t 



1 + Kn 



, fovj<£, 



= 0, iorj>£, 



1 + 



sfi + 4) 1 J'Cj'+tcj fV(i + «l)(i + «r 2 ) 



-(-e-i) 



y / (l + K ?)(l + K2 - 2 ) ' + , for j < £, 



= 0, for j > 



(9.3) 



Comparing with (|8.ip it is readily seen that these Kraus operators are linearly independent for exactly the same 
reason by which the quantum limited transpose channel was extremal. 
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IV. The composite V(k 2 ;0) oCi(«i;0), k 2 > 0, < ki < 1 

Kraus operators of this composite, which always corresponds to a noisy transpose channel (see Table IJ), are 
{T m (K 2 )B n (K 1 )}, < m,n < oo. We have 

oo 

T m (K 2 )B n (Ki) = ^2? mn \m- j)(n + j\, 

3=0 

x k{ ^1 - i^j , for j < m; 

= 0, for j > m. (9.4) 

It is immediately clear that T m (K 2 )B n (Ki) and T > (n2)B n ' (k\) are (trace-)orthogonal unless m + n = m + n . 

We may therefore divide the Kraus operators into orthogonal subsets determined by m + n = con- 
stant = N. Then linear independence will have to be established just within each subset Qjv = 
{To(k,2)Bn(ki), Ti(k2)Sjv-i(ki), • • ■ , Tn{k2)B (ki)}. It is seen from (|9.4j) that T (k2)Bn{ki) is a multiple of 
|0)(iV|, T 1 (k 2 )-Btv-i( k i) is a linear combination of |0)(-/V| and |l)(iV — 1|, T 2 (k2)B^_ 2 (ki) is a linear combina- 
tion of |0)(JV|, 1 1) (AT — 1|, and |2) (AT — 2|, and so on. Thus linear independence within f2jy follows as an immediate 
consequence of the fact that ^ nn ^0Vm,n. 



V. The composite Ci(k2;0)oD( Ki ;0), ki>0,0<k 2 <1 

This composite channel corresponds to a quantum-limited transpose channel (see Table [TJ) which we have already 
considered in much detail in Section III. The Kraus operators {B m (K 2 )T n (ni)} , < m, n < oo (which as a set should 
be equivalent to {Te(K 2 Ki)}, < t < 00), are 

n 

B m (K 2 )T n (Ki) = £L n \3 -m)(n- j\, 

j=m 




= 0, for 71 < to. (9.5) 

Thus B m (K 2 )T n (Ki) = if m > n. Further, all B rn {K2)T n {K\) with m = n correspond to a multiple of the vacuum 
projector |0)(0|, showing that these Kraus operators are not linearly independent. Even so, this is also a valid 
representation of the quantum- limited channel C\{k2', 0) o 2?(«i; 0) = T>(k2Ki', 0) . 



VI. The composite C2(«2;0) oD(/t 1 ;0), k% > 1, K\ > 



This composite channel corresponds, for all Ki,k 2 , to a noisy transpose channel, similar to the case of 2?(/t 2 ;0) o 
Ci(fti;0) considered earlier. The Kraus operators {A m (K 2 )T n (Ki)} , < to, n < 00 have the form 

n 

A m (K 2 )T n (Ki) = Y &nn\3 +m){n- j\, 

3=0 

Jm+jn.nn. fj-\ ( Ja J. ^ 



—j 



1 + K a 2 J , for j < n; 
0, for j > n. (9.6) 
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Linear independence of the Kraus operators can be established in a manner similar to the earlier case of V(k2',0) o 
Ci(ki;0). 

VII. The composite X>(k2; 0) o £2(^1; 0), K2 > 0, k\ > 1 

This composite is a quantum-limited transpose channel (see Table HJ, with Kraus operators {T m (K2)A n (Ki)}, < 
m,n < 00. Similar to the earlier quantum-limited case of Ci(ki;0) o 2?(k 2 ;0), these Kraus operators too are not 
linearly independent; they do represent the quantum-limited channel T>(k 2 ki]0), though. 

We conclude this Section with some further observations. It is seen from TableUthat to realise the noisy attenuator 
channel Ci(k; a), for the full parameter range < k < 1, a > 0, where a is a measure of the additional classical noise 
above the quantum limit, as the composite £2(^2; 0) Ci(ki;0), we have to solve K2«i = K and 2{k 2 — 1) = a for 
< /ti < 1, K2 > 1, and we have as solution k 2 = yl + a/2 > 1, «i = k/k 2 < k < 1. It should be appreciated that 
this realisation contains as special cases the noisy A\{a) for k = (i.e., Ki = 0), and the classical noise channel i? 2 ( a ) 
for k = 1. 

Similarly, the noisy amplifier C2{k;o), k > 1, a > can be realised through the same composite if we solve can 
K2K1 = k and 2k|(1 — k\) = a for < K\ < 1, K2 > 1. We have the solution «2 = \/ n 2 + a/2 > n > 1 and 
«i = k/k 2 < 1. Again the classical noise channel $2(0) is contained as the special case k = 1. 

Finally, to realise the noisy transpose channel P(/c; a), k > 0, a > as the composite T>(k2', 0) oCi(ki; 0), < Ki < 
1, K2 > 0, of quantum-limited channels we solve k — k 2 ki and a — 2k|(1 — k\) to obtain k 2 = \J n 2 + a/2 > k, Ki = 
k/«2 < 1- The same can also be realised as the composite £2(^2! 0) "D (ki;0), k 2 > 1, Ki > by solving k = K1K2 
and a = 2{k^ — 1), and we have K2 = yl + a/2 > 1, K\ = k/k 2 > 0. 

We may summarise some aspects of our consideration thus far in this Section in the following manner. 

Theorem 14 All the nonsingular noisy channels can be realised as the composition of a pair of quantum-limited 
channels. Equivalently, and as a consequence, each of the nonsingular Gaussian channel has an operator-sum repre- 
sentation in terms of a discrete set of linearly independent Kraus operators. 

The above assertion includes in particular the case of the classical noise channel #2(0). Further, the case A\(a) is not 
exempted from our consideration above, for it is just a special case, corresponding to k = of the noisy attenuator 
Ci(k;o). 

It is seen from Table U that the noisy singular case A2 (a) can be realised as a composite of two quantum- limited 
channels: either following or preceding the quantum-limited A2{0) by quantum limited Ci(-;0), C2(-;0), or 2?(-;0). 
Consequently, the Kraus operators will be indexed by one discrete and one real variable. Thus, the single quadrature 
noise channel i3i(a) is the only singular case that does not submit itself to our consideration above, in the sense that 
there seems to be no way of realising it as composite of a pair of quantum limited channels. 

While we have obtained in the Section Kraus representations for noisy channels with the aid of pairs of quantum- 
limited channels it is, of course, possible to obtain Kraus representation using unitary dilation of Section 2, with the 
ancilla in a thermal state rather than the vacuum state. But we believe our present approach has the advantage of 
leading to Kraus operators of extremely simple structure, in addition to the advantage of connecting the noisy case 
to the quantum-limited case. 

We conclude this section with three remarks, two of them are in respect of Table Q] while the third one is in the 
context of error correctability. 

Remark 1 : We have already noted in Sections III and V that the quantum- limited family T>(-) is self dual 
T>(k) ~ whereas the families Ci(-) and C 2 (-) are dual to one another : C\{n) ~ C 2 (k~ 1 ). The reader will 

recall that these duality relations are a consequence of the failure to be unital of these (trace-preserving) maps by 
just a multiplicative scalar. Remnants of these duality relations may be readily observed in Table HI The composite 
£>(•) oCi(-) is quantum-limited, the dual fact being that C2(-) oD(') is quantum-limited. The fact that C2O) °C2( - ) 
is quantum-limited is dual to the fact that Ci(-) o Ci(-) is quantum-limited. The fact that £>(•) o Ci(-) is a noisy 
conjugator has as its dual the fact that C2{ ) ° 2?(-) is a noisy conjugator. Finally, T>(-) o £>(•), Ci(-) o C2O, an d 
C2O Ci(-) arc self duals. Since the failure of the quantum limited A2 to be unital is nontrivial, A2 does not figure 
in any such duality relation. 

Remark 2 : In Table U we have considered the composition of pairs of quantum-limited channels under the 
assumption that the two channels are simultaneously in their canonical forms. In order to help the reader appreciate 
this remark it should first be emphasised that a typical composite of this kind, say 2?(«2;0) oCi(ki;0), should 
not stand for composition of two quantum-limited channels which are already in their respective canonical forms, 
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but rather to two channels picked one from either Sp(2,R) orbit or double coset. That is 2?(k 2 ;0) oCi(«i;0), 
for instance should stand for (([/(Si) o 2?(k 2 ;0) ° U(S2)) ° {U(Si) o Ci(ki;0) o (U(S 2 )), for arbitrary mctaplectic 
unitaries corresponding to Si, £2, S± , S 2 € Sp(2,R). The fact that two channels cannot in general be taken to their 
respective canonical forms simultaneously, i.e., with S ± = Si, S 2 = S2, brings out the nontriviality of the assumption 
underlying Table 1. When this assumption is lifted, TableUgets much enriched into Table 2 as shown in the Appendix. 

Remark 3 : Wc have noted in Sections III to VI that the Kraus operators {We} of quantum-limited Gaussian channels 
possess the property that the associated nonnegative operators w\Wt are simultaneously diagonal (in the Fock basis 
for £>(•), Ci(-), and C 2 (-) and in the position basis for A.-i). In the present Section we presented for each nonsingular 
noisy channel a discrete set of Kraus operators, say W mn , indexed by a pair of integer variables m, n and it can be 
readily verified in each case that the associated nonnegative operators W^ n W mn are simultaneously diagonal (in the 
Fock basis). For the noisy channel A2(a), it may be verified in the realization £>(•) o A2 and C 2 (-) o A2 [and not in 
_4 2 o £>(•), „4 2 °Ci(-), or A2°C2(-) ] that the relevant nonnegative operators are simultaneously diagonal in the position 
basis. Finally, the single quadrature classical noise channel Bi(a) being random unitary, the associated nonnegative 
operators are all multiples of unity. We may thus state 

Theorem 15 For every Gaussian channel it is possible to obtain a Kraus representation such that the nonnegative 
operators associated with the Kraus operators are all simultaneously diagonal. 

The above observation leads to the following remark on error correction. 

Remark on Error correction : The fact that W\ W/s are simultaneously diagonal in the Fock basis for all i for the 
channels D(k;0), Ci(k;0), C 2 (k;0) and their composites, imply in the view of the work of [39| that the Fock states 
could be used to reliably transmit classical information through this channel. In such a case the channel is viewed 
as a generalised measurement. In other words, any classical information encoded in Fock states and passed through 
these channels can be reliably retrieved by a restoring channel. 

X. CONCLUSION 

We have obtained operator-sum representations for all single-mode bosonic Gaussian channels presented in their 
respective canonical forms. Evidently, the operator-sum representation of a channel not in the caonical form follows 
by adjoining of appropriate unitary Gaussian evolutions before and after the channel. The Kraus operators were 
obtained from the matrix elements of the two-mode metaplectic unitary operator which effects the channel action 
on a single mode. The two-mode symplectic transformation in each case did not mix the position and momentum 
variables and this fact proved valuable for our study. The Kraus operators for the quantum-limited channels except 
the singular case were found to have a simple and sparse structure in the Fock basis. 

It was shown that the phase conjugation channels T>(k) and I?(k _1 ) are dual to one another, and the attenuator and 
the amplifier families Ci(k) and (^(k -1 ), k < 1 are mutually dual. The channels T>(n), Ci(k), and C 2 (k) were found 
to be almost unital; in the sense that the unit operator was taken to a scalar times the unit operator. The channel 
2?(1) was found to be bistochastic but not random unitary. The unitary bistochastic channels being the classical noise 
channels £>i and B2 and the trivial identity channel C\{1) = C 2 (l). 

In the case of the phase conjugation channel, the action in phase space was brought out explicitly through the 
action of the Kraus operators on the Fock basis. The attenuator channel resulted in the scaling of the diagonal weight 
function 4>(a) and the amplifier channel resulted in the scaling of the Husimi Q-function as expected. Further, the 
output of the channel with respect to classicality /nonclassicality was studied. It was found that the phase conjugation 
channel T>(k) and the singular channel A2 are classicality breaking while the attenuator channel Ci (n) and the amplifier 
channel C 2 (k) do not generate nonclassicality. 

The action of the channel in the Fock basis gave an insight into the fixed points of the channel. The action in the 
Fock basis together with the action of the channel in phase space led us to conclude that there is a unique thermal 
state which is an invariant state for T>(k), k < 1. Further it was shown that the vacuum state is the only invariant 
state for the attenuator channel Ci(k), and that there is no finite energy state that is invariant for either the amplifier 
channel C 2 (k) and the singular channel *4 2 . 

Using Choi's theorem, it was shown that all quantum-limited bosonic Gaussian channels are extremal. The Kraus 
operators of the phase conjugation channel was brought to a rank one form, thereby explicitly bringing out the entan- 
glement breaking nature of the phase conjugation channel. It was further shown that there is no finite rank operator in 
the support of the Kraus operators of either the amplifier or the attenuator channel, and this explicitly demonstrates 
that the quantum-limited attenuator and the amplifier families of channels are not entanglement breaking. The Kraus 
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operators of the singular channel Ai was also obtained in the rank one form thereby manifestly showing that this 
channel is entanglement breaking. 

It was shown that every noisy Gaussian channel (except the singular case Si (a)), can be obtained as the composition 
of a pair of quantum-limited channels as shown in Table [I] This in turn implies that apart from those compositions 
that involve Ai, there is a discrete set of Kraus operators for all the noisy Gaussian channels. Further, the nonnegative 
operators {W^Wi} were found to be simultaneously diagonal for all Gaussian channels. This throws light on the error 
correctability of these channels. 

In bringing out the semigroup structure of the amplifier and the attenuator families of quantum-limited channels, 
it was shown that interrupted evolution slows down both amplification and attenuation in a manner characteristic of 
the quantum Zeno effect. 
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Appendix : Composition of a General pair of Quantum Limited Channels 

Given two quantum-limited Gaussian channels whose (X, Y) matrices in the sense of Section 1 are respectively 
(Xi,Yi), (X 2 ,Y 2 ), with |detXj| = ft?, we have for the composite channel (X,Y) 

X = X1X2, Y = X2 Y\X 2 + Y2. 

If (X%,Yf), (X$,Y£) are the canonical forms of (Xi,Fi), (X 2 ,Y 2 ) in the sense of (TOj) . the most general (Jfi,Yi), 
(X 2 ,Y 2 ) should necessarily have the form (SiX°S 2 , SjY^S 2 ), {S 3 X$S 4 , SfY 2 °S 4 ), with S U S 2 ,S 3 ,S 4 € Sp(2,R), so 
that 

X = SiXiS 2 S3X 2 S4, 

Y = X® S3 S 2 YiS 2 SsX 2 S4 + S^Y® S4. 

_ T 

Our problem now is to classify the orbits or cosets under the unitary equivalence (X, Y) ~ (SiXS 2 , S2 YS 2 ), 
Si,S 2 € Sp(2,R). Basically we have to determine the determinants of X, Y in terms of the canonical parameters 
Ki,K2 of the constituent channels. While detX is independent of S±, S 2 , S3, S4, detY needs a careful consideration. 
It is this situation that Y > is the sum of two (positive) terms that breaks our analysis into two distinct cases. 

Let us first consider the case in which X® is nonsingular, so that both the terms of Y are nonsingular : this case 
obviously corresponds to the first twelve entries of Table HI1 With the choice S2 = S^ 1 the second term of Y becomes 
a multiple of the identity, and the first term of Y becomes a multiple of S3 Sj S2S3 € Sp{2, R), a positive symplectic 
matrix. We can now do a rotation S2 = R £ SO(2) C Sp(2, R) without affecting the second term in Y, so that 
the first term of Y becomes a multiple of diag(A, A -1 ), A, A -1 being eigenvalues of Sj S2 S2S3. DetY" can now be 
easily evaluated. Removal of the mandatory quantum- limited noise, as dictated by the value of detX [i.e., |1 — k 2 |, 
1 + ft 2 , or 1 depending on det X being positive, negative or zero, k 2 equalling det|X| ], then yields the classical noise 
indicated in Table fll] by the second argument of the composite channel. This procedure is the one used for the first 
twelve entries of Table HU 

The case of singular X2 is somewhat different. With the removal of S4 with the choice S 2 = S^ 1 , as in the earlier 

case, in the first term of Y the projection X® = ^ picks out the 1,1 element of the matrix S3 S 2 S 2 S3- With 

the eigenvalues of this matrix again denoted (A, A" 1 ) and assuming that is the rotation needed to diagonalizc it, 
the 1,1 element equals A cos 2 9 + A -1 sin 2 8 (a combination appearing in the last four entries of Table [II]). With this, 
the second term in Y is a multiple of identity whereas the first term is a multiple of A cos 2 8 + A -1 sin 2 8 times the 

projector ^ jM ; the determinant of Y can be readily computed, and removal of the quantum-limited noise gives the 

numerical value of the second argument of the composite channel. 

Finally, the two cases in the last entry of Table HI1 the case of the composite .42 (0) o A 2 (0), correspond to the 
two possible situations that may arise with X when both X® and X® are of rank one : X could cither be rank one 
or it could be a null matrix. Correspondingly, the composite ^2(0) ° -42(0) should be viewed as a member of the 
Ci(0; a) = -4i(a) family or the A 2 (a) family, the extra classical noise being 'the same' in both cases. 



36 



Ci(k 2 ;0) 


oCi(ki;0) 




Ci (k 2 ki; - nWi) 2 + - «?)(! - - A- 


-1)2 _ 




C 2 (k 2 ;0) 


oC 2 («i;0) 




C 2 (ac 2 ki; v/(k?«| - l) 2 + 4(k( - - 1)(A - A~ 


-1)2 _ 


- - 1)) 


^2\K2, 0) 


~ n Ct- . (\\ 


n 
Cl 


f K 2 Kl ; ^/ (2cc 2 — «f k| 


_1)2 +K 2 (1 _ K 2 )(K 2_ 1)(A _ 

for k 2 ki < 1; 


A J 


2 n .,2 2x \ 

— — K 2 KiJ 1 , 






C9 


f k 2 ki; a/(2k| — k^cc?. 


-1) 2 + K 2 (1- K 2 )( K 2 -1)(A- 
for k 2 ki > 1 


A- 1 ! 


2 - (k2k? - 1^ 

V n '2 ^1 A ^ i ) 


Ci(k 2 ; UJ 


„ /-i c t - . n\ 
o C 2 (ki; Uj 


/I 
Ci 


( k 2 ki; ^/ (fi; 2 ,^ 2 — 2k| 


+ 1)2 + K 2 (1 _ K 2 )(K 2_ 1)(A _ 

for K2K1 < 1; 


A J 


— (.1 - «2«lJ J, 






Co 

*-*2 


k 2 ki; ^/(kjkJ; — 2k| 


+ 1) 2 + K 2 (1-^ 2 )( K 2 -1)(A- 

for K2^1 > 1. 


A- 1 ! 


2 - f/eiK? - 1^ 
\<^2 ™1 L J J 7 


t>c • ni 


^ f>C ■ nl 




/ / { 9 9 i r» 9 
[ K2«l! V ( K 2 K 1 + ^ K 2 


1 1 \o i 9 / 1 i 9 \ / 1 1 9 \ / \ 
+ l) 2 + + Kf)(l + K|)(A - 

tor K2K1 <. 1; 


\— 11 


2 Ci t - 2 „2\\ 






c 2 


f K2K1', \/ (^2^1 ^^ 2 


4. 1)2 + K 2(l _|_ K 2)n + K 2 )(A - 
for K2^1 > 1. 


A- 1 ) 


2 - (^1^1 - 1)), 


£>(k 2 ;0) 


oCi(ki;0) 


V 


( v (2^1 ~~ K 2 K i 


+ 1) 2 + K 2 (1- K 2 )(1 + K 2 )(A- 


a- 1 ; 


2 -(1 + «!«?)) 


£>(k 2 ;0) 


oC 2 (aci;0) 




2? («2«i; 1/(1 + K 2 K i 


) 2 + ac 2 (^-1)(1 + k 2 )(A-A- 


1)2 _ 




C 2 (k 2 ;0) 


o£>(ki;0) 


V 


( k 2 kx; \J (2k| + «l K; i 


-1) 2 + K 2 (1 + K 2 )( K 2 -1)(A- 


a- 1 ; 


2_( 1 + K 2 K 2 ) j 


Ci(k 2 ;0) 


oO(aci;0) 




© ^« 2 ki; 1/(1 + ^2 K i 


) 2 +ac 2 (^ 2 + 1)(1- k 2 )(A-A- 


1)2 _ 


(! + «!«?)) 


Ci(/t 2 ) 


o^ 2 (0) 




A2 ( 


y/1 + — K 2 )(A — A -1 ) 2 — ] 


) 










■^2 (\/(2 


, 2 -i) 2 + ^ 2 ( K 2 -l)(A-A-i) 


2 1 






o^ 2 (0) 






, 2 + l) 2 + K 2 (l + K 2 )(A-A-i) 


2 - 1 


) 


A 2 (0) 


°Ci(m) 




-4 2 (\/l + (A cos 2 6» + A- 1 sin 2 <9)(1 - k 2 


-1) 




■r\2 \y ) 


u i^ 2 V^l J 




^(V 1 + ( A cos2 61 + A" 1 sin 2 6){kI - 1] 


— 1 1 




A 2 {0) 


oV(ki) 




A 2 (s/l + (A cos 2 9 + A- 1 sin 2 6») (1 + k 2 


"I) 




A 2 (Q) 


o^ 2 (0) 




A2 ( \/l + A cos 2 6> + A- 1 sin 2 6 - 1 j , 
Ci f 0; Vl + A cos 2 6» + A- 1 sin 2 6» - 1 


or 

)■ 





TABLE II: Showing composition of two quantum-limited Gaussian channels which are not necessarily in their respective 
canonical forms. It may be seen that the case A = 1 coresponds to Table |T] In the case of the first twelve entries the noise is 
always greater than what obtains in Table [I] But in the remaining four cases the noise can be either more or less, depending 
on the value of 8. 



We conclude with the following observations in respect of Table HI1 It is evident that in all cases A = 1 corresponds 
to the situation in which both constituents of the composite are already in their canonical forms, and the reader can 
verify that Table [TTI reduces to Table Q] in this case. And thus we see that A 7^ A -1 results, in the first sixteen entries 
of Table HH in classical noise which is always more in magnitude than the case A = 1. The last four entries of Table 
HI! are distinguished in this regard : since the classical noise depends also on the choice of 9, it can be either more or 
less than what obtains in the A = 1 case. 
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